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Abstract Liquid hold-up and gas phase pressure drop data are reported for mercury irrigation 
of a packed bed of | in. steel raschig rings with countercurrent flow of hydrogen and nitrogen. 


The static hold-up of mercury was found to decrease slightly with increasing liquid rates. 
This effect, not observed in aqueous or organic liquid irrigations, is attributed to an increased 
mobility of the liquid in the bed because of the greater energy associated with the flow of dense 
media such as a liquid metal. 


When nitrogen was the gas phose, the stacie hold-ups, dynamic hold-ups and irrigated liquid 
specific surfaces are all larger than those observed under comparable conditions with hydrogen 
as the gas phase. This is considered to be due to the greater interfacial tension of mercury and 
nitrogen, 


An empirical correlation is developed between the liquid hold-up, gas phase pressure drop, 
gas mass flow rate and gas density. It is suggested that this type of correlation would be of 
value for estimating the hold-up of molten metal or slag in high temperature systems, where direct 
measurement would be difficult, but where a room temperature analogue could be readily 
constructed. 


Résumé Des données concernant le hold-up dun liquide et la chute de pression dans la phase 
gazeuse sont relevées dans le cas de Virrigetion par du mercure d'un lit de garnissage d'anneaux 


de Raschig en acier df} in., avee un écoulement a contre courant @hydrogéne et dazote. 


Le hold-up statique du mercure décroit lentement avec les vitesses croissantes du liquide. 
Cet effet nest pas observe dans des irrigations aqueuses ou de liquides organiques ; il est attribué 
a une mobilitée accrue du liquide dans le lit, due a une énergie plus grande associce a Pécoulement 
du milieu dense tel celui dun metal liquide. 


Quand Pazote est la phase gazeuse les hold-ups statiques et dynamiques, les surfaces spécitiques 
irriguées par le liquide sont toutes plus grondes que celles observées observées avec Chydrogéne 
comme phase gazeuse, dans des conditions comparables. Ceci serait da a la plus grande tension in- 
terfaciale du mercure et de Vazote. 


Une corrélation empirique est développeée entre le hold-up liquide, la chute de pression dans 
la phase gazeuse, la vitesse de Péecoulement massique du gaz et la densité du gaz. Ll est suggéré 
que ce type de corrélation sera de valeur pour Vestimation du hold-up du métal fondu, ou dun 
laitier dans des systémes a haute température cl une mesure directe est difficile, mais of un 
analogue a la température ambiante peut étre rapidement construit. 


Zusammenfassung Der Fliissigkeitsinhalt und der gasseitige Druckabfall fiir die Berieselung 
eines Festbettes aus 6.5 mm-Raschigringen aus Stahl mit Quecksilber wird angegeben, wenn 


Wasserstoff oder Stickstoff entgegenstromen. 
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Der statische Flissigkeitsinhalt an Quecksilber fiel etwas ab mit steigendem Flissigkeitsstrom, 
wie es bisher bei wisserigen oder organischen Berieselungsflissigkeiten nicht beobachtet wurde. 
Die Ursache dieses Effekts ist dic gréssere Beweglichkeit der Flissigkeit in dem Bett wegen der 
hdéheren Energie, die mit der StrOmung einer Flissigkeit hoher Dichte, wie hier cines fliissigen 


Metalls, verbunden ist. 


Bei Stickstoff als strOmendem Gas waren der statische und der dynamische Fliissigkeitsinhalt 
und die beriesclte spezifische flissige Obertliche grosser als es unter vergleichbaren Bedingungen 
bei Wasserstoff beobachtet 
zwischen Quecksilber und Stickstoff gesehen. 


wurde. Die Ursache wird in der grésseren Grenzfliichenspannung 


Bezichung zwischen dem Flissigkeitsinhalt, dem gasseitigen 
Kine solche Beziehung kann fiir die 


Es wurde eine empirisehe 
Druckabfall, dem Gasstrom und der Gasdichte entwickelt,. 
\bsechitzung des Fliissigkeitsinhalts bei geschmolzenem Metall oder Schlacke in Hochtemperatur- 
systemen von Bedeutung sein, bei denen direkte Messungen schwierig, cin Analogieversuch 


bei Raumtemperatur aber leicht durchzufiihren ist. 


INTRODUCTION 


Mo.TeN METAL or slag irrigation of fixed beds 
under non-wetting conditions occurs frequently 
in pyrometallurgical operations. An analysis of 
such irrigations, for the purpose of evaluating 
heat and mass transfer rates, requires information 
on the volume of liquid in the bed, the liquid 
surface area and the condition of the liquid 
hold-up under operative conditions. 

To give an insight into the flow of molten slag 
over coke, countercurrent to a gas, GarpNer [1] 
measured pressure drops and hold-ups for water 
irrigating a packed coke bed under non-wetting 
conditions and indicated that the data were to 
enable interpretation of future heat and mass 
transfer studies on the molten slag system. 

A similar motive prompted this investigation, 
as information was required to analyse mass 
transfer data obtained when molten lead was the 
irrigating liquid phase and nitrogen containing 
zine vapour at an elevated temperature was the 
countercurrent gas phase. 

To simulate the flow of molten lead over the 
packing of an absorption column, mercury was 
used as the liquid phase with hydrogen and 
nitrogen as the countercurrent flowing gas phase. 

Although these studies were primarily aimed 
at evaluating the dynamics of molten lead flow 
over a packed bed, by studing the mercury 
irrigation of the same packing material, it is 
that the 
application to other liquid metal or heavy media 


considered discussions have general 


irrigations. 


THEORETICAL ASPECTS 


(i) Fluid flow through fixed beds 


Single phase fluid flow through packed columns 
has been reviewed by Excun [2]. 

The equation developed by Ercun [2] on the 
basis of numerous previous reported data for 
fluid flow in granular beds was 

AP g (1 € U% 
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i.e. the pressure losses are caused by simultaneous 
kinetic and viscous energy losses. The first term 
on the R.HLS. of equation (1), representing the 
energy the most important 
factor at low flow rates and of little significance 


Viscous loss, is 
at high flow rates, whereas the second term in 
the right-hand side of equation (1), representing 
the kinetic energy losses, is the most significant 
at high flow rates and of little importance at very 
low flow rates. As the gas flows involved in this 
investigation were fairly high, the viscous energy 
losses could be assumed negligible compared to 
kinetic energy losses, and the form of equation 
(1) simplified ; 

AP 


L &. a’ p | a 


Excun [2] assigned a value of 1-75 
{( ee) D, | for the term a’, for the single 
fluid phase flow through a packed bed. Under 
the conditions of these studies, the bed porosity 
« altered according to the hold-up of liquid in the 
packing as also did the size parameter D,. The 
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term a’ in this investigation was therefore con- 
sidered to be some function of the total hold-up 
of the liquid in the packing, and the true velocity 
(V) was used rather than the superficial velocity 
(U) in the modified equation (3). 


AP 
a (hy) 


x’ (pV) 


The method originally proposed by BLakr 
[3] and used subsequently by Carman [4] for 
correlating pressure drop data in packed beds, 
by an approach analogous to the established 
procedure of plotting the friction factor as a 
function of Reynolds number for fluid friction in 
pipe systems, was modified by Morr [5] to allow 
a more direct comparison to be made between 
porous media and pipe flow. The forms proposed 
by Morr are presented in equations (4), (5), 
(6), (7). 

1U p 
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As suggested by Garpner [1] the types of 


correlation, employing the surface area per unit 
packed volume, are suitable for examination of 
low when the packing is irrigated, as they allow 
the total surface area exposed to the gas flow to 
be estimated from a measurement of the gas phase 
pressure drop. The unknown S can be found by 
cover plotting the data on a previously determined 
correlation in which the dry surface area is known, 
Such an evaluation of the irrigated total specific 
surface would be useful for establishing the correct 
gas phase Reynolds number, but in itself it does 
not indicate the contribution to the specific 
surface, which is due to the liquid phase, as the 
specific surface obtained by such a_ procedure 
represents both solid and liquid surface presented 
to the gas flow. 


(ii) Liquid hold-up in packed beds 

The void space in the packing, which is occupied 
by liquid under operating conditions, must be 
taken into consideration when the bed porosity 
is determined. The total volume of liquid in the 
packed bed under operating conditions is known 
as the total hold-up. The total hold-up (h,) is 
considered to be the sum of the static hold-up 
(hs) and the dynamic hold-up (hp). The static 
hold-up is defined as the liquid in the packing 
which does not drain from the packing when the 
liquid supply to the column is discontinued. The 
dynamic hold-up represents the liquid which will 
drain from the packing when the liquid supply 
is cut off. The relation between the three hold-ups 
is given by equation (8). 


hp = hg + Np (8) 


Several investigators have measured or 
estimated liquid hold-ups. ‘Two extensive studies 
Eicixn and Wess [6] and 


Jesser and Excin [7] but there was no complete 


were reported by 


analysis of the three different types of hold-up 
for packed columns, until the recent work of 
SuvuLMan et al. [8] and Garpnexr [1]. 

These investigators postulated that the total 
One 
part of the hold-up was considered to be due to 


hold-up may be divided into two parts. 


accumulated semi-stagnant liquid, residing at 
places which have a natural retentive effect 
upon the liquid, such as points of contact between 
the packing elements. This liquid was considered 
to be semi-stagnant but not entirely at rest, the 
flow being so slow that its volume was governed 
by other effects such as the surface tension of the 
liquid. Suutman et al. [8] considered this com- 
ponent of the hold-up to be equal to the static 
hold-up. GARDNER [1] recognized the existence of 
further terms contributing to the semi-stagnant 
component. A small quantity was added to the 
static hold-up to allow for liquid passing in and out 
of the sites of static hold-up and a_ further 
allowance was made for the pressure drop across 
the 
liquid, 

The flowing part of the hold-up was considered 
the dynamic Garpner [1] 


developed an empirical expression allowing the 


bed creating new sites of semi-stagnant 


to be hold-up. 
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flowing part of the hold-up to be analysed in a 
logical form. 

GARDNER [1] and Sutman ef al. [8} both 
reported the static hold-up as being independent 
of the gas and liquid rate. The static hold-up 
was found to be a function of the physical proper- 
ties of the packing surface. 

In the past, it had been assumed that the 
dynamic hold-up was independent of gas rate up 
to the flooding point and so dynamic hold-ups in 
some of earlier studies were determined at zero 
gas rate. However, both Suutman ef al. [8] and 
Garpnek [1] showed that this was not the case, 


as the dynamic hold-up was found to increase with 


increase in gas rate, even before the loading point 
was reached. In the non-wetting case, it was 
suggested that the effect of gas rate on the 
dynamic hold-up was that the liquid was retarded 
by frictional effects or that it was dispersed into 
a larger number of thinner trickles. 

Besides the bed 
evaluated, hold-up studies are also important in 
the interpretation of mass transfer data, because 
of the different significance, which can be attached 
to the components of the total liquid hold-up, 


allowing porosity to be 


in so far as their being effective, in the mass 
transfer process. This 
was discussed in detail by Suutman ef al. [9], 
when the effective interfacial distinct 
from the total interfacial area, was analysed. 


very important aspect 


areca, as 


EXPERIMENTAL 


Apparatus and procedure 


A small glass column was packed to a height of 


1-75 ft with } in. steel raschig rings. The volume 
of rings in the column was obtained by weighing 
the rings prior to their addition to the column. 
The volume of the packed height of the empty 
glass column was determined by addition of a 
tared quantity of water, The relevant dimensions 
of the column are given in Table 1 and a line 
diagram of the experimental apparatus is shown 
in Fig. 1. 

The mercury feed to the column was from a head 
tank through a pipe provided with both a plug 
cock and a needle valve. The needle valve was 
used for flow control and the plug cock served 


as the quick acting valve required for measuring 
the mereury hold-up. The mercury flow rate was 


Table 1. 


Dimensions of packed column 





1-75 in. 
0-01665 ft? 
1-75 ft 

231 mi 

S26 ml 

O72 
(595-hy) 826 
Packing material } in. } in. 
steel raschiy rings 


Column diameter 

Free cross-sectional area 
Packed height 

Volume of rings 

Volume of empty bed 
Dry bed porosity 
Irrigated bed porosity 
20 gauge 





H 


A 


a 


a 


™ 


ee | 

1 

; 
ae, o 


Fic. 1, Diagram of apparatus. 

A-Stoekcock on inlet liquid. B-Stopeock on outlet liquid. 
C-Liquid run-off tube. D-Storage bowl. E-Raschig 
ring packing. F-Mercury head tank. G-—Variable lute, 
H- Receiving vessel. J-Manometer. K- Mercury entrain- 
ment trap. L-Needle control valve. M-Venturi meter. 
N-Rotameter. P- Exhaust to atmos- 
phere. Q-Compressed gas connexion. 


O- Drainage cock. 
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Table 2. Pressure drop and hold-up for mercury irrigation 





Mercury Static Dynamic Gas Total Bed Gauge True gas Pressure 
rate hold-up hold-up hold-up hold-up porosity pressure velocity drop 
(Ib ‘min) hs hyp he hy € at base of Ve across bed 
(ml) (ml) (ml) (ml) column (ft sec) AP 
(em Hq) (em H,O) 
Nitrogen gas phase 

1-50 112-9 0: 31° 153-9 
18-00 101-2 < 226-2 
516 108-6 140-6 
10-00 100-2 ; : 165-2 
18-84 92-9 175-9 
a3B4 103-2 : 134-7 
2-41 105-5 ° 130-0 
10-56 97-2 77° 184-7 
14-41 96-2 WO! M7 
44 103-2 i- 164-7 
3-06 103-2 22:s oY 139-2 
6-89 93-9 j : 134-9 
12-56 100-2 { 202-2 
2-59 98-2 i 2 126-2 
3-22 101-2 of 162-2 
135-75 96-2 5 154-2 
3-28 105°: i 1215 
5-16 102°: 123-2 
16-88 85°. ; 146-5 
84 Ot : 121-9 
3-62 : 110-2 
14°31 ; 52-5 143-3 
4-72 Os: 22 120-2 
2-03 , 110-2 


Hydrogen gas phase 
15-56 5 177-5 
8-06 f 131-5 
17-06 5 151-3 
384 “{ : 112-9 
sn) . f 125-9 
14-84 ° 139-8 
74 119-9 
6-038 i 117-2 
18-34 3 BHO 
10-84 125°8 
5-25 . 110-3 


- +) = 


uw st st 


20-47 -! 146-5 
3-41 101-8 
8-97 116-5 

18-78 5 | 156-5 


— ® 


8-09 118-9 
6-69 115-9 
13-59 135-8 
6-28 2 113-8 
2-78 111-2 
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obtained by weighing the discharge from the 
column. 

Nitrogen or hydrogen were metered to the 
column by a venturi meter for the higher flow 
The 


gas was admitted at the bottom of the column 


rates and a rotameter for the lower flows. 


with manometer tappings located above and 
below the packed bed. The pressure drop across 
the bed was measured using a simple manometer. 

The base of the glass column was provided with 
a bowl of suflicient storage capacity to collect 
the mercury hold-up when the supply to the 
column was cut off, 

The experimental technique was as follows: 


(a) When the 
desired gas and mercury flow rates were 


steady state conditions at 
reached, the variable lute was adjusted 
so that the level of the liquid in the 
run-off tube was at the calibration mark 
and the pressure drop across the bed was 
recorded, 
Stockcocks A 
simultaneously, but the gas supply to the 
column was left unaltered. 


and B were then closed 


The liquid which drained from the packing, 
the dynamic hold-up, was determined by 
opening stopcock B and draining off the 
the the 
liquid level in the run-off tube was again 


mercury into receiver H, until 
at the calibration mark. 
The stopeock B was closed once this level 


had 


to the column was then turned off. At 


been reached and the gas supply 
very high gas flow rates, a further quantity 
of mereury drained from the column, which 
the same as the 
This 
the total hold-up was termed the “ gas 


hold-up.” 


was measured in way 


dynamic hold-up. contribution to 


The static hold-up was determined by 
weighing the column at the conclusion of 
each run and subtracting the weight of the 
column before irrigation. 


Experimental results 


Using the foregoing procedures, the pressure 
drops and hold-ups for mercury irrigation of the 


packing with hydrogen and nitrogen flow were 
determined and the results presented in Table 2. 
The mercury hold-ups in the absence of gas flow 
but with hydrogen and nitrogen atmospheres 
are given in Table 3. 

Pressure drops were also measured across a 
previously mercury irrigated bed and across the 
dry packing before irrigation and are shown in 
Tables 4 and 5 respectively. 

As all data were obtained under approxi- 
mately the same ambient conditions, a constant 
temperature of 20°C was used throughout in 
evaluating the physical properties of the gas 
phases. The true gas velocity takes into account 
the actual void space in the bed under the operat- 
ing conditions and the density evaluated at the 
mean pressure in the bed. 


DISCUSSION 
(i) Static hold-up 


Previous investigators, who have reported 
static hold up data in aqueous and organic liquid 
systems, have observed the static hold-up as 
being constant for the particular liquid and 
packing. However, the results of this investigation 
show that the liquid metal irrigation rate has a 
slight but definite effect on the static hold-up, in 
that the static hold-up decreases with increased 
liquid rate. This effect can be seen in Fig. 2. 
Although the liquid rate effects the static hold-up, 
no consistent effect could be attributed to the gas 


flow rate although the results at zero gas flow 


rate do appear slightly higher than the other data. 

The other interesting feature of Fig. 2 is that 
the static hold-up for nitrogen is displaced above 
the hydrogen data. This effect is thought to be 
due to a difference in interfacial tension between 


the mercury and the two gases. Perry [9] gives 
the interfacial tension of mercury and nitrogen 
as higher than the interfacial tension of mercury 
and hydrogen, 496 dynes/em and 
$70 dynes /cm at 15°C and 19°C respectively. The 
increase of static hold-up with an increase in 
surface tension is consistent with the finding of 
Suu MAN et al. [10], who studied a surface tension 
range of 23 to 86 dynes ‘cm by the addition of 
wetting agents to water. 

A probable explanation of the observed effect 


being 
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Fie. 2. Mercury irrigation statie hold-up. 


Nitrogen atmosphere, no flow. 
@ Nitrogen flow. 


which the liquid rate had on the static hold-up 
in this investigation is that liquid metals or other 
density liquids possess much greater kinetic 
energies than water when irrigating a packing, 
and therefore the possibility of dispersing or 
“ knocking-on ” semi-stagnant pockets of liquid 
is much greater and this effect increases with 
increase liquid rate. The significance of this would 
be that the semi-stagnant pockets of liquid, 
observed in packings irrigated with water, would 
be very much more mobile in the case of a dense 
liquid irrigating the packing. Under these 
conditions, the whole concept of effective inter- 
facial area in mass transfer would have to be 
reviewed, as the increased mobility of the dense 
liquids could well result in the static hold-up 
being almost as effective as the recognized flowing 
component, the dynamic hold-up. 


(ii) Dynamic hold-up 

The dynamic hold-up for mercury irrigating 
the packing is not purely a function of liquid 
rate, and therefore supports the recent results 
of GarpNer [1] and Suvutman et al. [8], who 
showed that the earlier practice of assuming 


mt 


Dynomic hold-up, 


Hydrogen atmosphere, no flow. 


A Hydrogen flow. 














8 2 16 
Mercury rote, ib /min 


Mercury irrigation dynamic hold-up. 


Nitrogen absence of gas hold-up 
@ Nitrogen presence of gas hold-up. 

Hydrogen absence of gas hold-up. 
A Hydrogen presence of gas hold-up, 





N. A. WaRNER 


the dynamic hold-up to be solely related to the 
liquid rate was incorrect. 


Table 3. Mercury hold-up in absence of gas flow 





Total 
hold-up 
(ml) hp hy 

(ml) 


Static hold-up Dynamic 


Mere 
lercury rate he hold-up 


(ib min) 


(ml) 


Hydrogen atmosphere 
20-49 SO-8 a” 
11-41 2-4 30 
5-69 OG 16 
ooo Os 13 
Nitrogen almosphere 
sz re 
1O6-5 pt) 
100-6 va. 


M20 LRU 





Table 4. Pressure drop across previously mercury- 
irrigated column 
ooo 


Previous irrigation hy = 98 mi and « 





True gas Pressure drop 
across bed AP 


fom TO) 


Gauge pressure 
at base of 


column (om Hyg) 


velocity 
Py (ft sec) 
Hydrogen gas phase 
83 
1 


cm & 


os «ss 
as - & 


— 


Nitrogen gas phase 
RB 
30 
oo 
70 
a8 





Fig. 3 is a plot of dynamic hold-up versus the 
mereury irrigation rate. The dynamic hold-up 
is clearly dependent on other factors besides the 
liquid rate. Increased gas rates, in general, 
increased the dynamic hold-up. The reasons for 
The high 


as rates correspond to the points in Fig. 3 in 


this were outlined by Garpner [1]. 
g 
which there was a “* 
hold-up plotted is the liquid which drained from 


gas hold-up.” The dynamic 


the packing before the gas flow was reduced and 


therefore excludes any “ gas hold-up” contribu- 


tion. 


Table 5. Pressure drop across column before 


irrigation € ‘72 





Pre assure drop 
across bed AP 
(om HO) 


Gauge pressure True gas 


at base of velocity 
column (om Hye) I (ft sec) 


Nitrogen gas phase 


Las phase 





(iii) Reynolds number—friction factor correlation 
The experimental data of this investigation 
analysed the 
GARDNER [1] to see if the total specific surface 
could be readily obtained, thus allowing the 


were according to method of 


correct Reynolds number to be assigned to the 


gas phase. 
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As the specific surface appears in the denomina- 
tor of both the Reynolds number and the friction 
factor, the most convenient method of evaluating 
the specific surface is to assign a value of unity 
to the dry packing, thus enabling the irrigated 
specific surface to be calculated as a ratio of the 
dry surface of the packing, which may be calcu- 
lated, if the loss of area at the points of contact is 
neglected. The calculation is essentially a trial 
and error process of substituting values of the 
ratio of irrigated specific surface to the dry 
specific surface until the irrigated experimental 
point cover plots the curve representing the dry 
data, which is used as the datum for all caleula- 
tions, 

The curve representing the dry data, in which 
the specific surface is taken to be unity, is 
obtained by substituting 8 
(7) and (5) to obtain 


1 into equations 


AP 2g « WVepe 
L pV fa 
as the functions to be plotted. 


The but 
previously irrigated packing, were also calculated 


data for zero mercury flow rate, 


in the same form, as it was expected that the two 
sets of data would only vary by a constant ratio 











and would, therefore, plot parallel to each other 


on a log-log graph. The logarithmic plot of 


2 9a 4 
ars < versus at Ve hes 
L V,? Pp be 


is shown in Fig. 4 in which the solid line is fitted 
to the experimental results. 

The unexpected feature of Fig. 4 is that the 
specific surface of the dry and previously irrigated 
packing is the same. For both cases a single line 
represents all the data plotted. The conclusion 
may be reached for this particular instance, where 
the previous irrigation was such that a_ total 
hold-up of 98 ml remained in the packing, that the 
increase in surface area due to the presence of 
the liquid hold-up is exactly compensated by the 
decrease in the surface area of the dry packing 
exposed to gas flow, 

The full line in Fig. 4 was then used as the 
datum line for evaluation of the irrigated specific 
surfaces. Values of 

AP 2g _— 4 pe 

VipL ie Sa 

eP t 
were calculated for all the irrigated data, and 
the results plotted on a log-log graph, in which 











ig. 4. 
Nitrogen before irrigation. 
Hydrogen before irrigation, 


%Ppse 4 
- 


Dry and previously irrigated packing. 


@ Nitrogen zero Hg rate 
A Hydrogen zero Hg rate 


























O Nitrogen Hg irrigated. 


the full line of Fig. 4 is drawn as the datum line 
as shown in Fig. 5. 

Fig. 5 displays an unusual characteristic in 
that the hydrogen data, except for the highest 
effects 


were operative, exhibit specific surface of the 


Reynolds numbers, where “ gas hold-up” 
same order or perhaps slightly less than the 
whereas the nitrogen data suggests 


dry packing, 
that the irrigated specific surfaces are quite 
significantly higher than the dry packing specific 
surface. 

However, the the 


gas rate investigated, which are not plotted in 


data for lowest nitrogen 


Fig. 5, showed a reversal in trend, as the irrigated 


specific surfaces were less than the dry packing 


specific surface. As the pressure drops were 
measured with a simple manometer the error 
involved in the measurement of the pressure drops 
at low gas rates was such that the low gas rate 
data were of doubtful accuracy. 
decided not to plot these points but to record 
that they showed a reversal of trend. To be con- 


sistent all other pressure drops less than 1-5 em 


Hence, it was 


H,O were also not plotted. 
A possible explanation of the significantly higher 
specific surfaces of the nitrogen data, even at gas 


rates where there was no “ gas hold-up,” is that 
differences in interfacial tension are the cause of 


the observed behaviour. Thus it would appear 


Vepte io" 


ff 


Mercury irrigated data. 


@ Hydrogen Hg irrigated. 


that, when nitrogen is the gas phase, the static 
hold-up and the specific surface for mercury 
irrigation are greater than the corresponding 
values for hydrogen as the gas phase, because of 
the higher interfacial tension of nitrogen with 
mereury. 

In cases where there is a “gas hold-up,” 
the specific surface would be expected to be 
somewhat greater than the other data. The group 
of points in Fig. 5 for values of 


4Vepe x lo* 


a 


greater than 15 for nitrogen are cases in which 
there were appreciable gas hold-ups. 

The two broken lines drawn parallel to the 
datum line in Fig. 5 represent specific surfaces 
30 per cent above and below the dry packing 
specific surface. The only data which show 
specific surfaces significantly outside these limits 
are the nitrogen data, in which the gas hold-up 
has increased the specific surface. 

(iv) The and 


pressure drop 


relation between total hold-up 


The type of correlation suggested by equation 
(3) appeared to be satisfactory for analysing 
the data, 
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\+3 5x10 he 


P 























Fic. 6. 

Nitrogen Hg irrigated. 

A Hydrogen Hg irrigated 
@ Nitrogen zero Hg rate. 


AP 
a (hr) 


AP. [x (hy) ]. for a given value of 


r’ (V2 p) (3) 


AP 


AP. ”" 


x (hy) = 


if x (hy) 


a + bhp, where a, b, ¢ are constants. 


Since AP = AP., when hy = 0 


a (hy) 1 + bhy (10) 


The constants b and e¢ were determined by 
fitting an equation to the curve obtained when 
selected values of x (h7), calculated from equation 
(9), were plotted against hy. Using this procedure 
the form of a(hy) given in equation (11) was 
obtained. 


a (hy) 35 x 10° hf) (11) 


All the data of this investigation could be 


empirically correlated by a straight line on a 


log-log plot of 


AP 


(1 +35 = W°hF 


_ +0 
Kp 


Empirical correlation of all data. 


A Hydrogen zero Hg rate. 
@ Nitrogen before irrigation. 
Iiyvdrogen before irrigation. 


versus V2 p as shown in Fig. 6. 

The correlation shown in Fig. 6 was developed 
for mercury irrigation, dry packing and previously 
mercury-irrigated packing. However, it seemed 
reasonable to expect the same correlation to be 
applicable to other liquid phases besides mercury. 
therefore undertaken the same 
packing with water irrigation at various liquid 
and gas flow rates. The packing was then treated 
with paraflin wax to make it non-wetting to 
water and similar non-wetting 
conditions were made. 

The correlation developed from the mercury 
studies applied reasonably well to the water 
irrigated data, and hence the proposed correlation 
appears to be satisfactory, at least as a first 
approximation, for estimating the liquid hold-up 
on the same packing for any liquid system, by 
measuring the gas phase pressure drop, the mass 
flow rate and the density. 

The other important feature of the correlation 
developed is that the state of the irrigating 
liquid, whether at rest as in the previously 
irrigated data or under actual irrigation, appears 
to make little difference to its applicability. 


Tests were on 


tests under 
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The correlation also holds equally well when the 
column is loaded, as the data in which there were 
appreciable the bed are 
represented by the same line as obtained for the 
column below the loading point. This is considered 


“gas hold-ups in 


to be indicative of the close inter-relationship of 


pressure drop and hold-up, as most other forms 
of analysis show anomalies above the loading 


point. 


Estimation of liquid hold-ups 


Although Suvutman et al. [10] have investi- 
gated the properties of the liquid phase, which 
affect the hold-up of aqueous and organic systems, 
the extrapolation of their results to a liquid 
metal of vastly different physical properties 
would be unwise. 

More work of a fundamental the 
flow of liquid metals or similar dense media is 


nature on 


required before a_ reliable assessment can be 


made of the exact effect which the physical 
properties of the liquid phase have on the liquid 


hold-up. 

In the absence of 
suggested that the empirical method just outlined 
could be used for estimating liquid hold-ups in 
beds when the experimental difficulty precludes 
direct measurement, as for example, in high 
temperature liquid metal or molten slag irrigations, 
If the empirical relationship between the total 
liquid hold-up and the gas phase pressure drop 
can be developed on a room temperature model 
of the system, the extension of the correlation 
to the high temperature system would enable 
the total hold-up at zero irrigation rate (i.e. the 
static hold-up) or the total hold-up with irrigation 


such information, it is 


to be estimated from the observed gas phase 
pressure drop across the bed, the gas temperature 
and = the The dynamic 
hold-up under the particular conditions would 
then be the difference of these estimates. 


gas mass flow rates, 
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NOMENCLATURE 


characteristic dimension of the particle ft 
dimensionless 
ft see® 


Fanning friction factor 
acceleration due to gravity 
gravitational conversion factor 
dynamic hold-up mil 
gas hold-up mi 
static hold-up mil 
total hold-up mi 
length of packed bed ft 
pressure drop across packed bed, Ib ft? in 
equations (1, 2, 6, 7), but for convenience 
in the dimensional equation (3) and the 
em H,O 
pressure drop across the dry packed bed 
em H,O 


dimensionless 


empirical correlation of Fig. 6 


Reynolds number 
surface area per unit packed volume of 
bed ft? ft® 
superficial fluid velocity ft/sec 
true fluid velocity ft sec 
true gas velocity ft 
fractional voids or bed porosity 
fluid density 

fluid viscosity 

function of total liquid hold-up 
function of gas inertia. 


Ib ft4 
Ib ft sec 
x (hy) 
x vi? p) 


REFERENCES 


[1] Garpner G. C. Chem. Engng. Sci. 1956 5 101. 
[2] Eroun 8S. Chem, Engng. Progr. 1952 48 89. 


[3] Brake F.C. Trans. Amer. Inst. Chem. Engrs. 1922 14 415. 

[4] Carman P. C. Trans. Inst. Chem. Engrs. 1987 15 150. 

[5] Morr R. A. Some Aspects of Fluid Flow pp. 242-256. Edward Arnold, London 1950, 
[6] Ever J.C. and Wetss F. B. Industr. Engng. Chem, 1989 31 435, 


[7] Jesser B. W. and Exverw J.C. Trans. Amer. Inst. Chem. Engrs, 1943 34 2 
[8] Suunman H. L., Untiarcn C.F. and Weiis N. Amer. Inst. Chem. Engrs. J. 1955 


277. 


247. 


[9] Perry J. H. Chemical Engineers Handbook (3rd Edition). p. 363. MeGraw-Hill, New York 1950, 


[10] 


Suutman H. L., Uneaicn C. FL, Weis N. and Proutx A. Z. Amer. Inst. Chem. Engrs. J. 1955 


259. 


160 





Chemical Engineering Science, 1959, Vol. 11, pp. 161 to 182. Pergamon Press Ltd., London. Printed in Great Britain 


The absorption of zinc vapour in molten lead 


Countercurrent absorption in a packed column 


N. A. Warner 


School of Metallurgy, The University of New South Wales, P.O. Box 1, Kensington, N.S.W., Australia 
(Received 25 August 1958: in revised form 11 May 1959 


Abstract Experiments are described in which nitrogen containing zine vapour is counter- 
currently contacted with molten lead in a small column packed with | in. mild steel raschig rings. 
The heights of overall gas phase transfer units for the absorption of zine vapour in molten lead 
using this packing are reported for a range of gas and lead temperatures and flow rates. 

It is suggested that the scrubbing of gases containing zine vapour, as in the recently developed 
zine blast furnace, should be considered as an absorption not a condensation. It it demonstrated 
in these studies that zine vapour can be continucusly removed from a gas phase by absorption 
in molten lead under conditions such that both the gas and lead temperatures are well above 
the zine dew point of the gas mixture, 

\ fundamental equation, for analysing the absorption of zine vapour in molten lead in 
terms of conventional gas absorption theories is derived using the additive diffusional resistance 
concepts. Room temperature analogues of the zine absorption in molten lead packed column 
enable evaluation of the flow characteristics of the irrigating liquid lead and the mass transfer 
properties of the packing under non-wetting conditions. 

For the experimental conditions of these studies, the absorption of zine vapour in molten 
lead is shown to be controlled by the diffusional resistance in the gas phase and can therefore 
be compared with other gas phase controlled absorptions such as ammonia in water from a dilute 
ammonia—air mixture. To correlate zine and ammonia absorption data, a new quantity, the 
effective hold-up ratio is introduced. The results of the zine absorption studies suggest that 
the total liquid hold-up is effective in this process because of the mobility of the static hold-up, 
when a liquid metal irrigates a packed bed. 

The successful correlation of both the ammonia absorption in water and the zine absorption 
in lead data indicates that the absorption of a metallic vapour in a molten metal at an elevated 
temperature can be analysed according to conventional gas absorption theories and also 
shows the value of room temperature analogues for experimentally difficult high temperature 


mass transfer processes. 


Résumé—L auteur décrit des expériences dans lesquelles de lazote contenant de la vapeur 
de zine est mis en contact en contre-courant avee du plomb fondu, dans une petite colonne a 
yarnissage contenant des anneaux raschig en acier doux d°} in, Les hauteurs d‘unité; de transfert 
de la phase gazeuse globale pour labsorption de la vapeur de zinc sur du plomb fondu en utilisant 
ce garnissage sont données pour un intervalle des températures du gaz et du plomb ct des vitesses 
d'écoulement. 

L/auteur suggére que lépuration de gaz contenant de la vapeur de zine soit comme dans 
les hauts fourneaux a zine récemment développés, considérée comme une absorption et non 
comme une condensation. Cette étude montre que la vapeur de zine peut étre continuellement 
éliminée dune phase gaz par absorption par du plomb fondu dans des conditions telles que les 


températures aussi bien du gaz que du plomb soient bien au-dessus du point de rosée du zine du 


mélange gazeux. 

Une équation fondamentale est définie pour Panalyse de labsorption de la vapeur de zinc 
par le plomb fondu en fonction des théories habituelles de labsor, tion des gas, en utilisant les 
concepts supplémentaires de la résistance a la diffusion. Des analogues a une température 
ambiante pour labsorption du zine dans une colonne a garnissage par du plomb fondu permettent 
d'évaluer les caractéristiques découlement du plomb liquide irrigant et les propriétés de transfert 


de masse du garnissage dans des conditions de non-mouillage. 
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Pour les conditions expérimentales, absorption de vapeur de zine par du zine fondu est 
réglée par la résistance 4 la diffusion dans la phase gaz et peut done étre comparée avec d'autres 
absorptions réglées par la phase gaz, telles que absorption de lammoniac par Peau a partir d'un 
mélange dilué air-ammoniac. L introduction dune nouvelle quantité, le rapport de hold-up 
(rétention) effectif permet de coordonner les données d'absorption du zine et de lammoniac. 
Les résultats des études sur labsorption du zine suggérent que le hold-up liquide total est effectif 
dans ce procédé & cause de la mobilité du hold-up statique quand un métal a état liquide irrigue 
un lit de garnissage. 

La concordance des données sur labsorption de Tammoniac par leau et sur labsorption 
du zine par le plomb montre que labsorption d'une vapeur métallique a température élevée, 
peut étre analysée d’aprés les theories ordinaires sur labsorption des gaz, et donne la valeur 
des analogues & température ambiante pour les transferts de masse a température élevée dont 
lexpérimentation est difficile. 


Zusammenfassung Es werden Versuche beschrieben, in denen Zinkdampf enthaltender 
Stickstoff mit geschmolzenem Blei im Gegenstrom in Beriihrung gebacht wird in einer kleinen 
Kolonne, die mit 6,5 mm-Raschigringen aus Fluss-Stahl gefiillt ist. Die Gesamt-Cbertragungs- 
hdhen der Gasphase fiir die Absorption von Zinkdampf in geschmolzenem Blei in dieser Packung 
werden fiir einen Bereich der Gasund Bleitemperaturen und Stromungsmengen mitgeteilt. 

Es wird vorgeschlagen, die Reinigung des Zinkdampf enthaltenden Gases, wie in dem kiirzlich 
entwickelten Zinkhochofen, als Absorption und nicht als Kondensation anzusehen, Die Unter- 
suchungen zeigen, dass Zinkdampf kontinuierlich aus der Gasphase durch Absorption in ge- 
schmolzenem Blei entfernt werden kann, wenn Gas-und Bleitemperaturen deutlich tiber dem 
Taupunkt des Zinks liegen. 

Unter Benutzung des Begriffs der additiven Diffusionswiderstiinde wird fiir die Absorption 
des Zinkdampfes in geschmolzenem Blei eine Grundgleichung in Ausdriicken der gewOhnlichen 
Absorptionstheorie abgeleitet. Analogieversuche bei Raumtemperatur in der gleichen Kolonne 
erlauben die Auswertung der Fliessgréssen des herabrieselnden fliissigen Bleis und der Stoff- 
iibertragungseigenschaften der Packung bei Nichtbenetzung. 

Unter den experimentellen Bedingungen dieser Untersuchungen zeigt sich die Absorption 
von Zinkdampf in geschmolzenem Blei bestimmt durch den Diffusionswiderstand in’ der 
Gasphase. Sie kann daher etwa mit der Absorption von Ammoniak in Wasser aus einer verdiinnten 
\mmoniak-Luft-Mischung verglichen werden. Zur Korrelation der Zink-und der Ammoniak- 
absorptionswerte wird cine neue Grosse, das Verhiltnis der effektiven Fliissigkeitsinhalte, 
eingefiihrt. Die Ergebnisse der Untersuchungen tiber die Zinkabsorption zeigen, dass der gesamte 
Flissigkeitsinhalt wirksam ist wegen der Beweglichkeit des statischen Flissigkeitsinhalts, wenn 
ein fliissiges Metall ein Festbett berieselt. 

Die erfolgreiche Korrelation der Ammoniakabsorption in Wasser und der Zinkalsorption 
in Blei zeigt, dass die gewOhnliche Absorptionstheorie auch fiir die Absorption eines Metall- 
dampfes in geschmolzenem Metall bei erhéhten Temperaturen giiltig ist; auch der Wert von 
Analogieversuchen) bei Raumtemperatur§ fiir Stoffiibertragungsprozesse bei cxperimentell 
schwierigen hohen Temperaturen wird nachgewiesen. 


INTRODUCTION molten lead. Continuous removal of both heat 
Tue development of the zine blast furnace by and zine from the furnace gases is obtained by 
the Imperial Smelting Corporation at Avonmouth, pumping molten lead through a scrubber and an 
England, has been recently outlined by Morcan external cooling and phase separating vessel, 
{1}. The success of the new process depends from which the liquid zine produced is continu- 
primarily on the recovery of zinc from a gas phase ously withdrawn and the lead returned to the 


containing a relatively low concentration of zine scrubber. 

vapour and a considerable amount of carbon The data of Table 1 were given by MorGan [1] 
dioxide. To prevent re-oxidation of the zine for the countercurrent scrubbing of zine blast 
vapour on cooling, the hot gases are shock- furnace gases with cool liquid lead saturated with 


chilled by direct contacting with showers of zine. 
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Table 1. Typical countercurrent scrubber operating 
conditions (based on figures by MorGan [1)) 





Gas phase 
5°, Zn at 1,000 C to be cooled to 450°C 


Lead phase 
Circulation rate 
Inlet temp. 

Inlet composition 
Outlet temp. 
Outlet composition 


$20 ton /ton Zn produced 
wc 

2-02°,, Zn (Saturated) 
570 ¢ 

2-26°,, Zn 


Saturation at outlet temp. +-90°,, Zn 





MorGan [1] and the zine blast furnace patent 
literature refer to the gas—liquid lead contacting 
stages as condensers. It is proposed that such a 
nomenclature the fundamental 
mechanism involved, as does the statement made 
by MorGan [1] when describing these units: 


overlooks 


* Rapid cooling tends to promote fog forma- 
tion but the seriousness of this aspect is 
reduced by (a) the high superheat of the 
gases, which have a dew point of about 
650°C and leave the furnace at 900-1000 C : 
and (b) the presence of a rain of unsaturated 
lead in which droplets of zine are formed, 
but are themselves scrubbed by the intense 
rain of lead.” 


This implies that discrete drops of zine are 
formed in the contacting stage, because the gas 
temperature in the vicinity of the lead falls below 
the dew point and it is these droplets of zinc 
which are collected by the scrubbing action of the 
lead showers. 

The transfer of zine from the gas phase into the 
liquid lead is clearly a diffusion controlled mass 
transfer operation. The zinc must diffuse from the 
bulk of the furnace gases through the laminar 
gas film on the gas-liquid lead interface and must 
then diffuse into the bulk of the lead phase. 
Hence at no stage would the zine be present in 
particulate form. 

Zine droplets would be formed in the bulk of 
the gas phase, only if the bulk gas temperature 
fell below the dew point of its zinc vapour content. 
As the dew point of the inlet gases is approximately 
650°C, the gas inlet temperature 1,000°C and the 
outlet lead temperature 570°C, it is very diflicult 


Countercurrent absorption in a packed column 


to conceive the bulk gas temperature as falling 
below the dew point, when the temperature 
drops associated with the transfer of heat from the 
hot gases to the cooler lead are taken into account 
and when it is realised that the zinc content of the 
gases will decrease approximately exponentially 
* condenser,” with the 


. 


as the gas traverses the 
corresponding decrease in the dew point tempera- 
ture of the gas mixture. Local cooling and zinc 
mist formation at the lead surfaces would also be 
unlikely, because in the gas film surrounding the 
lead surfaces a large zinc partial pressure gradient 
would be expected at high rates of mass transfer. 

A more appropriate nomenclature for the gas 
liquid lead contacting stages would be to call these 
units “ absorbers,” as it would appear that the 
fundamental unit operation involved is basically 
gas absorption. 

The aim of the present investigation was to 
study the mechanism of the absorption of zinc 
vapour in molten lead and to see if the conven- 
tional gas absorption theories, developed for 
non-metallic systems at or near room temperature, 
could be applied to the absorption of a metallic 
vapour in a liquid metal at an elevated tempera- 
ture. 

Metallic vapour-liquid metals systems have not 
previously been analysed from this point of view 
and thus there have been no similar investigations 
reported in the literature. 


EXPERIMENTAL 
When the 


initiated, the laboratory disk column, as introduced 


experimental programme was 
by Sreruens and Morris [2], was considered to 
be the most desirable apparatus for studying the 


continuous countercurrent absorption of zinc 


vapour in molten lead. Because of difficulties 
associated with the flow of molten lead over the 


A small 


packed column, whilst not possessing the advant- 


disks, this column proved unsuitable. 


age of known interfacial area, was finally chosen. 

The guiding considerations in the design of the 
experimental apparatus were that the absorption 
would have to be carried out on a relatively small 
that of the 
apparatus for periods up to one day would be 


scale and continuous operation 


required to collect the necessary data. 
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Although the absorption column itself was only 
of a small size, the ancillary systems, essential 


to its operation, presented numerous design 


problems because of the high temperature of 


operation. 

Basically the experimental apparatus consisted 
of a small packed column irrigated with molten 
lead and up which a stream of nitrogen and zine 
vapour could be passed with provision for measur- 
ing the gas and lead phase temperatures, flow 
The 
components and layout of the apparatus are shown 
in Fig. 1. 


rates and zine compositions. essential 


Description of apparatus 

(i) Liquid lead system. The molten lead supply 
to the column was from a large insulated tank 
containing a small submerged centrifugal pump 
and needle control valve and provided with an 
inert nitrogen atmosphere. Intermittent checks 
on the lead discharge rate into a casting wheel 
located at the bottom of the absorption column 
The lead 
gin. O.D. 


18 gauge stainless steel tube and was heated to 


enabled control of the lead flow rate. 


discharged from the pump into a 


the desired temperature as it went through this 
tube by passing a heavy electrical current through 
the pipe 
350A were required to heat the lead to the column 


and its contents. Approximately 
operating temperature. The heavy current source 
1000 A, 12-5 V Selenium Rectitier 


The passage of heavy current 


used was a 
Unit. 


through the tube also served to preheat the 


Plating 


tube to above the melting point of lead, before the 
lead pump was switched on. 

Direct heating by the passage of heavy current 
provided a very convenient method of keeping 
other components at the required temperature, 
The superheater and gas analysis sampling lines 
in the gas circulation system were kept well above 
the dew point of the zinc vapour by this method, 
The low voltage, high current electrical circuit 
used is shown in Fig. 2. 

(ii) Gas circulation system. Nitrogen was used 
The 


circuit operated on a closed system with two 


as carrier gas for the zine vapour. 


gas 


sliding vane rotary blowers in parallel providing 
A wet 


the necessary boosting for recirculation. 


gasholder filled with light oil was connected to 
the suction side of the boosters to maintain a 
positive pressure of 2 in. w.g. on the suction side 
of the system. The infiltration of air into the 
system through any leaks was thus prevented. 
Any leaks in the system were made up by supply 
of nitrogen from the gasholder. The nitrogen 
supply to the gasholder was by an automatic 
solenoid operated valve from a compressed 
nitrogen source. 

(a) Preheater. 


temperature of approximately 600°C, a 


To preheat the nitrogen to a 
heater 
was constructed consisting essentially of a tube 
furnace inside which were six passes each 2 ft 6 in, 
long of 2 in. O.D. stainless steel tube. 

(b) Carbonizer. The commercial grade nitro- 
gen make-up to the system was contaminated 
with oxygen and therefore a unit was designed, 
so that the oxygen content of the nitrogen was 
carbon monoxide before 


converted to passage 


to the zine boiler. This was achieved by inter- 
posing, in the line between the preheater and the 
zine boiler, a tube packed with carbon granules 
maintained at approximately 1,000 C, This unit 
was designated the carbonizer. 

As the nitrogen was recycled, the problem of 
oxygen contamination would have been eli- 
minated, if there were no leaks in the gas circula- 
the 


removal of oxygen would have been substantially 


tion system requiring make-up, because 
complete by gaseous oxidation of the zine after 
one pass through the zine boiler, 

(c) Zine boiler. Preheated nitrogen was passed 
over a crucible containing liquid zine at or near 
the boiling point (906°C at latm). This unit 
was termed the zine boiler. The power input to 
the boiler determined the rate of vaporization to 
give the desired zine partial pressure in the 
nitrogen stream. 
the service 
required and the extreme corrosive nature of 


Because of high temperature 
liquid zine at these temperatures, external heating 
of a metal boiler containing a crucible of zine 
was precluded, The back pressure on the boiler 
was approximately 5 Ib in® gauge and therefore 
the boiler had to be contained in a metal vessel 
to prevent excessive leakage and to facilitate 
connexion to the other metallic components of 
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Fic. 2. Low voltage, heavy current electrical circuit. 
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the gas circulation system. The solution to these 
considerations was to build the furnace inside a 
mild steel vessel and to insulate the space between 
the furnace and the metal walls so that the metal 
wall temperature was below 600°C with the clectric 
resistance furnace in the zine containing atmos- 
phere and under the pressure of the unit. The 
contructional details of the zine boiler are shown 
in Fig. 3. 

(d) Superheater. To ensure that there was no 
condensation of zinc in the gas line between the 


zinc boiler and the absorption column, a super- 
heating section was installed to bring the gas 
temperature to 700-800°C before it entered the 
absorption column, For the mass flow rates 
required, the necessary superheating was obtained 
in a four foot length of jin. O.D. stainless steel 
tube at 850°C set in between the boiler and the 
absorption column. <A current of approximately 
250 A was passed through this tube, to maintain 
the superheater tube at the required temperature. 

(iii) The packed absorption column. The 
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Fic. 3. 

A-Morgan sillimanite crucible N3907. 
D-Alumina cement. 
gas manifold, J-Silica tube. 
N-Insulating brick. 
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B-M organ sillimanite crucible N1938. 
EK. Zine charge. F—Power input insulating tubes. 
K-Inlet gas sampling line. 
O-Vermiculite insulation. 


The zine boiler. 


C—* Kanthal ” resistance winding. 
G-—Heavy current connexion. H_-Outlet 
L.-Superheater. M-Gas inlet from carbonizer. 
P-Steel wool. 





The absorption of zinc vapour in molten lead. Countercurrent absorption in a packed column 


principal dimensions of the packed column are Table 2. Details of the packed column 
summarized in Table 2. 





: Internal diameter 1-74 in. 
The shell of the absorption column was 


fabricated from a 2 in. O.D. 10 gauge stainless 
steel tube. Details of the shell construction are gauge steel raschig rings 
shown in Fig. 4. The tower was packed to a Packed height 1-75 ft 


5 
height of 1-75 ft with jin. mild steel raschig Bed porostty O73 
(before irrigation) 


Empty cross-sectional area 0-01654 ft® 


Packing material } in. jin. x 20 


rings. The raschig rings were machined by 
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Fic. 4. The packed absorption column, 
A-Packed bed } in. MUS. raschig rings. B-Column shell. C—Column furnace. D-Packing support. E-Lead tempera- 
ture probe. F- Lead collecting reservoir. G-Lead inlet nozzle. H-Lead inlet from pump. J-Gas inlet temperature 
probe. K Gas outlet temperature probe. L-Gas inlet superheater. M-Inlet gas distributor. © N-—Liquid seal. 
O-Gas outlet from column. P-Outlet gas sampling line. 


167 





N. A. WaRNER 


parting off jin. lengths of | in, O.D. 20 gauge 
mild steel tube. 
the column was recorded to allow calculation of 
the bed porosity. 

Lead and gas temperature probe tubes were 
welded to the bottom section of the shell. Through 
these probe tubes, stainless steel thermowells 


The weight of rings added to 


made from 3-5 mm O.D. stainless steel tube were 
inserted into the lead and gas streams. The lead 
located in a small reservoir 
placed directly under the 
Portion of the lead flowing down the column 
collected in this 
temperature to be measured. The gas temperature 


was located directly beneath the 


thermowell was 


packing support. 


reservoir, thus enabling its 
thermowell 
reservoir and therefore contact with lead flowing 
in the column was precluded. 

A manometer pressure tapping was located 
beside the gas temperature probe at the bottom 
of this column. This pressure tapping was also 
shielded from the lead flow by placing it directly 
underneath the reservoir used for measuring the 
lead temperature. 

No connexions were made to the top section 
of the column shell of a 
2kW (240 V) nichrome resistance furnace on a 
the shell to 
loss from the 


to allow installation 


column 
the 


silica tube around 


the 


3 in. 
balance heat walls of 
column. 

The top gas temperature probe and the mano- 
meter the 


sealing ring welded to the top cover plate of the 


pressure tapping were located in 
column after the furnace was placed in position, 
The top lead temperature was measured by a 
thermowell located in the liquid lead feed nozzle. 

To the top cover plate of the column were 
welded the lead inlet tube, the gas outlet and 
combined primary zine condenser and entrained 
lead trap and the outlet gas analysis sampling 
line. 

The bottom cover plate of the column contained 
the gas inlet distributor. The gas distribution to 
the bottom of the column was through 5 » | in. 
diameter ports. 

To prevent zine 
heating was provided around the gas distributor at 
the bottom of the column by a coal gas burner 
formed in two semicircular sections so that the 


condensation, additional 


circumference of the distributor was evenly 
heated. A gas burner was also provided on the 
top cover plate. 

The lead discharge pipe from the column was 
submerged in a reservoir of molten lead in which 
an overflow was positioned to allow continuous 
removal of the from the The 
liquid seal thus formed was capable of holding 
The 


into 


lead column. 


column pressures up to 51b in® gauge. 
overtlow from the liquid seal was cast 
ingots, 

To prevent liquid lead from flowing back into 
the gas circulation system in the case of a blockage 
in the lead discharge from the column, a relay 
device was installed in the absorption column to 
interlock lead level in the column to the centrifugal 
pump. Once the lead inside the column built up to 
a dangerous level, the lead pump automatically 


off. A 32V 


absorption column. 


probe was located inside the 
Once 
between the insulated probe and the wall of the 


cut 


contact was made 


liquid lead, a change-over 
electromagnetic relay de-energized the coil of an 
cut off the 


supply of current to the motor driving the pump. 


column by type 


electromagnetic contactor, which 


The contactor was also arranged with thermal 
overload protection for the motor. 


Analytical procedures 
The 


nitrogen containing zine vapour was analysed for 


(a) Gas phase zine analysis system, 
its zine content before and after leaving the 
absorption column. Samples of gas were drawn 
off from the inlet and outlet gas streams by piping 
sample lines to the suction side of the gas boosters, 
The sample of gas for analysis on the inlet 
side of the absorption column was taken from the 
outlet manifold of the zine boiler. The sample of 
gas for outlet zine analysis was taken directly 
from the top cover plate of the absorption column. 
Both inlet and outlet sample lines were of 3 in. 
O.D. stainless steel tube connected electrically 
together in series. The lines formed the circuit 
for the passage of the heavy current (approxi- 
mately 250A) used for keeping them above 
the dew point of the zinc vapour. Temperature 
control of the lines was by insertion of the 
required resistance in series with the circuit, 
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The sample lines were connected by “ ermeto ” 
couplings to removeable stainless steel headers 
inside resistance furnaces also maintained well 
above the dew point of the zinc vapour. Connected 
to the headers and protruding out of the analysis 
furnaces were welded Zin. O.D. stainless steel 
tubes, each containing (at the cold end) a small 
glass wool filter unit for collecting the condensed 
zine formed, as the gas cooled inside these tubes. 
These filter units were connected by detachable 
j in. O.D. copper tubes to a manifold provided 
with a plug cock for each tube, allowing a sample 
to be drawn through the desired analysis tube. 
Interposed between the manifolds and the blower 
were 30 ft® hr oil filled gas meters on both the 
inlet and outlet analysis systems. Orifice meters 
were located on the outlet of the gas meters to 
enable the rate of sampling during a run to be 
observed and controlled at the desired rate. 

Provision was made in the construction of the 
analysis headers for installing seven tubes in both 
the inlet and outlet analysis headers. The inlet 
outlet identical. 
Stainless steel all welded construction was used 


and analysis headers were 
throughout, to enable the tubes to be leached 
with nitric acid and to prevent undue oxidation at 
the high temperatures used. At the conclusion 
of the experimental run on the apparatus, the 
headers were removed from the furnaces, the zine 
leached out of each tube individually with hot 
5 N HNO, and determined volumetrically using 
sodium versenate. 

(b) Chemical analysis. After an initial separa- 
tion of lead by a two stage precipitation with 
sulphuric acid, the zine in the lead effluent from 
the column was determined volumetrically using 
the disodium salt of ethylenediamine tetracetic 
acid (sodium versenate) under buffered conditions 
at pH 10 with Eriochrome Black T as the indicator. 
Complete separation of lead and zine was not 
required initially as the small amount of lead 
remaining in solution was titrated with sodium 
versenate with the zinc ions complexed by cyanide 
addition. Once the lead end point was reached, 
the addition of formaldehyde destroyed the zinc 
complex rendering the zine available for titration 
with sodium versenate. 

The method outlined was checked with synthetic 
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mixtures of lead nitrate containing small known 
amounts of zine and was found reliable. 

A primary separation of lead and zine was not 
required in the gas phase zine analysis. Because 
of the low lead partial pressures, the zine collected 
in the analysis tubes was only slightly con- 
taminated with lead and hence it was possible to 
titrate this lead directly, whilst the zine was in 
the cyanide complex, before releasing the zine for 
titration. 


Operating procedure 

The gas circulation system was purged with the 
blowers on, until the residual oxygen in the system 
was less than 0-5 per cent, before the power to 
the various furnaces in the circuit was turned on. 
Complete oxygen at this stage was considered 
unnecessary, as also was the use of special oxygen- 
free nitrogen. As the nitrogen in the system was 
recycled, the final removal of oxygen from the 
system relied on the reaction with the hot charcoal 
bed in the carbonizer and also on reaction of any 
CO, or O, with the zine vapour formed later 
when the zine boiler had reached operating 
temperature. The gas in the system, under the 
final operating conditions, was therefore considered 
to be oxygen free, as the only make-up to the 
system was to allow for leaks and this small 
make-up passed through the carbonizer before 
entering the zine boiler and the absorption 
column. 

When the gas circulation system had reached 
steady state conditions for the column tempera- 
tures required, the lead pump was switched on 
and the flow adjusted by the submerged needle 
valve. After 10 to 30 minutes of pumping, the 


column temperatures became fairly steady and the 


run was commenced. (Minimum number of phase 
throughputs before run, liquid lead = 25, 
gas = 500). For a period of 10 to 30 min, depend- 
ing on the lead flow rate, the following run 
procedure was followed : 


(i) All the lead was collected and weighed. 
(ii) The total gas flow over the test period 
was read from the gas meter. 

The inlet and outlet analysis samples 
through the appropriate 


(iii) 


were drawn 
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analysis tubes, tagged for the particular 


run and the sample volumes noted. 

The rate of sample withdrawal was noted 
by observing the orifice meters on both 
the inlet and outlet analysis systems. 
The pressure at the base of the column 
was read. For the lower gas through- 


puts, it was sometimes necessary to 
increase the column pressure by slightly 
closing a valve on the gas line before the 
meter, to allow analysis samples to be 
drawn off at the desired rate. 

(vi) The rotameter reading and pressure at 

the blower were noted. 

(vii) The pressure drop across the packed 

column was measured. 

(viii) Samples of lead for analysis were obtained 

by directly collecting the discharge from 

the liquid seal in a steel crucible approxi- 

mately T}in. high fin. diameter. 

Usually six samples were obtained evenly 

spread over the test period. 

The total time of the test period was 

noted the the run 


marked on the chart of the sixteen point 


and duration of 
temperature recorder which recorded all 
relevant operating temperatures. 


Although the gas analysis system was designed 
to allow seven runs on the apparatus before 
removal of the analysis headers for zine analysis, 
it was found that only three runs could be made 
during a one day campaign. The usual procedure 
was to allow the gas system to cool down over- 
night, with the nitrogen still connected to the 
gasholder, after the completion of three runs, and 
to repeat the procedure the following day for the 
remaining three runs, before dismantling the gas 
analysis headers for analysis. The remaining gas 
analysis tube in the inlet and outlet assemblies, 
not used in the series of runs, served as a blank in 
the zine analysis. 


Experimental results 


The experimental data obtained using the 
foregoing procedures are tabulated in Table 3. 
Runs 4B, 4C and 4D were discarded because of a 


leak discovered at the conclusion of the experi- 


ments in the inlet analysis tube No. 4. This leak 
was not present in Run 4A, as the analysis tubes 
used in the A series of runs were replaced at the 
conclusion of this series by modified forms to 
enable easier dismantling. Run 1D was discarded 
because the automatic cut-off on the lead pump 
was actuated during this run and therefore a 
complete series of data was not obtained, 

For series A, the gas inlet zine mole fraction 
was calculated from the experimentally deter- 
mined gas outlet and lead analyses, assuming a 
zine balance of 100 per cent. For Series B, C 
and D, the gas inlet mole fraction was obtained by 
chemical analysis. The reason for this difference 
was that in Series A, the first set of runs on the 
the procedure 
different to the later runs in that inlet gas was 


apparatus, experimental was 
not analysed during the test period but measured 
directly after the test. This procedure was later 
considered unwise, as operating conditions tended 
to alter and therefore the carly procedure did 
not give a true indication of conditions during the 
test period, 

The outlet lead was analysed for zine in Series 
A and B. This enabled a check to be made on the 
accuracy of the analytical procedures including 
the efliciency of the zine traps and the gas and 
lead 
balances for these runs. A zine balance could not 
Runs 1A and 2A 


initial experimental difliculties in these two runs 


flow systems, by calculating zine mass 


be established on because 


prevented samples of the inlet gas being taken. 
A and B 


These results clearly 


The zinc mass balances for Series 
Table 3. 
indicate the reliability of the experimental data, 
particularly when the experimental difliculties 


are shown in 


associated with the investigation are considered, 
As a result of the excellent zine mass balances 
obtained, the time consuming zine analysis of the 
outlet lead was not undertaken for Series C and D. 
In these series, the inlet zine concentration in the 
lead and the inlet and outlet gas phase zinc 
compositions were determined experimentally. 
The zinc in the outlet lead was calculated assum- 
ing a zinc mass balance of 100 per cent. Although 
the outlet lead was analysed in Series B, the values 
tabulated in Tables 3 were calculated 
assuming 100 per cent zine balance, as the gas 


also 
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phase analyses were considered more reliable 
than the determinations of zine in the outlet lead. 


Evaluation of over-all transfer units 


By definition, for dilute systems * the number 
of over-all gas phase transfer units” (Nog) is 
given by equation (1): 

¥y 
dy 
‘ (1) 
Ye 

For a particular absorption column of packed 
height Z, the number of transfer units and the 
height of a transfer unit are related : 


Z Hog Noe 


If the equilibrium line is linear, Nye can be 
calculated from the terminal concentrations : 


" Ye 


Nog 
(Y — Ye). 


(3) 
and if y, = 0 


Nog = In (2 (4) 
Y2 
For the small column used in this investigation 
experimental difficulty precluded direct measure- 
ment of the terms required to rigorously evaluate 
equation (1); only the terminal values were 
the 
evaluated using two assumptions. 


measured and intermediate conditions 

An exponential decrease of the gas phase zine 
concentration from the bottom to the top of the 
column was assumed, As the equilibrium zine 
partial pressures for the alloys formed in these 
studies were usually small, the assumption made 
was unlikely to introduce any serious error. 

The second assumption made was that the 
lead temperature was linearly related to the 
position in the column. As the exit lead tempera- 
ture was in no case greater than 1-35 times the 
inlet temperature, the more straightforward linear 


relationship was used rather than an exponential. 


Under the conditions of this investigation, only 
a slight error was introduced if Nog was evaluated 
by equation (3), using the terminal conditions 
only. As the difference between the two methods 
of evaluating Nog was small, the simpler procedure 
was used throughout to evaluate Np, and hence 


Hog. The values of Nog and Hog obtained using 
equation (3) are given in Table 4. The equilibrium 
zine partial pressures used in the calculations were 
based on the free energy change equations of 
Keiiey [3] and Lumspen [4], respectively, for 
the vapour pressure of pure zinc and the zine 
activity coeflicient in lead—zine alloys. 


Table 4. Calculated values of Nog and Hog 





Series A 

Run BA WA 
iz | a oa 
Ono O41 


Nog 
Hog (tt) 


Series B 
Run 2B 3B 
256 370 


O70 O48 


Nog 
Hog (ft) 


Series C 

Run 2c 3c 
1-58 1-78 
113 10) 


Nog 
Hog (tt) 


Series D 

Run 2D 
2:58 
0-68 


Nog 





Appitive DirrustonaL Resistance ANALYSIS 


For the countercurrent absorption of zine 
vapour in molten lead in a packed column, the 
temperature of the liquid lead changes appreciably 
as it flows through the absorption column, and 
the equilibrium conditions alter according to the 
position in the column because the zinc activity 
coeflicient and the vapour pressure of pure zinc 
vary greatly with temperature. 

As the absorptions of this investigation were 
carried out in very dilute alloys, Henry’s law 
could be assumed to apply for each temperature 
level within the absorption column, Hence the 
equilibrium zine partial pressure over the alloys 
formed by absorption, for a particular temperature, 
are linearly related to the zine concentration in 
the liquid lead alloy; i.e. the mole fraction (y,) of 
zine in the gas phase, which is in equilibrium with 
a zine mole fraction of (x) the liquid lead phase 
can be given by the following equations : 
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Pan © Yan Pan 


PZn YZn 2 (6) 


* 2 P, ) 


The equilibrium conditions in a system in which 
molten lead is contacted with zinc vapour can there- 
fore be represented, for the case of very dilute alloys, 


by a series of straight lines of slope (pz, yz,)/ Pr 


If the equilibrium 
this light, the 
fundamental equation (7) can be derived for the 


for each lead temperature. 
conditions are considered in 
absorption of zinc vapour in molten lead, using 
the additive diffusional resistance concepts. 
"1 uv 
dy _ G,, dy 
I Katy—y) Pjkaty 
Ys Wy 
2 *, 
4. G | (p Zn Y%n) dy 
ad d Py kya Pm (y 
v 


G,, 


“== 


Y,) 


yy 


This equation can be simplified somewhat if the 
individual film coeflicients are assumed constant 
throughout the column (strictly only correct for 
isothermal dilute 
rewritten in terms of transfer units as shown in 


systems) and the equation 


equation (8). 
Fa « 
G,, i, | ‘Pa Yan) dy 


Pry y “ 


Hy 


Ln N OG « Ye 


¥9 

As the value of Hyg is proportional to the 
over-all resistance to mass transfer, equation (8) 
can be considered as a statement of the addition 
of diffusional resistances of the individual phases. 
The first term in the R.HL.S. of equation (8) is a 
measure of the resistance to mass transfer in the 
gas phase and the second term accounts for 
resistance in the liquid lead phase. 

Interpretation of the experimental results in 
terms of equation (8) was not possible by further 
experimentation on the high temperature system 
and hence two room-temperature analogues of 
the lead—zine absorption column were constructed. 
The extension of the data derived from the low 
temperature models to the absorption of zinc 
vapour in molten lead required a critical examina- 
tion of the methods of correlating gas and liquid 
film mass transfer data. The empirical correlation 
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methods are therefore briefly reviewed in the 
appendix, 


Room temperature absorption model 


To develop correlations for estimating the 
height of a gas film transfer unit and the height 
of a liquid film transfer unit, ammonia and carbon 
dioxide were absorbed in water using the same 
tower packing as used in the zine absorption 
So that in. this 
temperature model would simulate absorption of 


experiments, conditions low 
zine vapour by molten lead, the tower walls and 
packing were made non-wetting, with respect to 
the irrigating liquid, by treatment with paraffin 
wax. 








ft at 25°C 


Nn 
oO 


True 4, 





- 
+ 














0 


Liquid rote, ib moles /hr 


hic. 5. Carbon dioxide absorption in water. True heights 


of liquid film transfer units, 
@ Non-wetting 21 in. bed. 
© Non-wetting 10-5 in. bed. 


The experimental techniques used in the room 
temperature absorption model were the same as 
used by Sterwens and Morris [2] and Tay.or 
and Rowers [5] for calibrating a disk column 
As these and similar 
techniques have been used frequently by other 


for physical absorption. 


investigators, the experimental methods con- 
tained no novel features and therefore the end 
results only are presented in Figs. 5 and 6. The 
heights of transfer units plotted have been 
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Gos rote, ib moles / hr 


Ammonia absorption in water. True heights of 
gas film transfer units. 

1-32 Ib mole hr) 

bed. 

bed. 


(Constant water rate 
@ Non-wetting 
oO Non-wetting 


21 in. 
10-5 in. 








Dynamic hold - up ratio 
The effect 


Ammonia absorption in water. 

dynamic hold-up ratio. 

(Constant air rate = 0-0705 Ib mole hr) 
© Non-wetting 10-5 in. bed. 


A Wetting 10-5 in. bed. 


corrected for end effects. The same correction of 


+8 in. of additional packed height was required 


for both absorptions to bring into line data 
derived from beds of 10-5 and 21 in. packed height. 

The effect of liquid rate on the height of a gas 
film transfer unit is given in Fig. 7, which shows 
that H, varies with H, °”. The choice of the 
effective hold-up ratio for studying liquid rate 
effects in gas phase mass transfer is discussed in 
the appendix. For the ammonia absorptions the 
effective hold-up ratio is considered to be the 
dynamic hold-up. 


Estimation of molten lead hold-up and velocity 


A second room-temperature analogue of the 
zinc absorption column, employing mercury as 
irrigating liquid metal with hydrogen and nitro- 
gen as the countercurrent flowing gas phase, was 
used to estimate the hold-up and average liquid 
velocities. The total hold-up of liquid lead under 
the operating conditions for each zine absorption 
run was estimated from the observed pressure 
drop across the column, the gas mass flow rate 
and the mean gas temperature and pressure in 
the bed. The estimation was a trial and error 


process of selecting values of total hold-up so 
that the experimental data cross-plotted the 


empirical correlation developed previously by the 
author [6] for mercury irrigation of the same bed. 
The estimated total liquid lead hold-ups for all 
runs are given in Table 3. 

As only the total hold-up of lead could be 
estimated using this procedure, the dynamic 
liquid lead hold-up was determined as_ the 
difference of the estimated total hold-up under 
the particular operating conditions and _ the 
estimated total hold-up in the absence of irrigation. 

The static hold-up of liquid lead in the packing 
was obtained by observing the pressure drop 
across the column, the gas mass flow rate and the 
mean gas temperature and pressure inside the bed 
immediately after the completion of an experi- 
mental run with the lead pump shut off. 
Examination of the value of hold-up thus 
obtained revealed that, prior to Run 5A, the 
hold-up of liquid lead was below “ saturation 
level’ and it was only after the high liquid lead 
rate of Run 5A that the static hold-up sites in 
the packing were completely occupied by semi- 
stagnant liquid lead. After Run 4A, the estimated 
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static hold-ups remained fairly constant at 
approximately 180 ml of liquid lead. 

The average linear velocity of water and 
mercury flowing under non-wetting conditions 
in the packing was experimentally determined by 
recording the time taken for the liquid phase to 
flow measured 
‘ saturated ” packing. 


over a height of previously 


As the liquid velocity was required to estimate 
the gas-liquid relative velocity in the packed bed, 
the high gas velocities associated with the zine 
absorption data enabled the lead velocity to be 
adequately estimated for the purpose of this 
investigation by assuming the lead velocity to be 
equal to the mercury velocity at the same volu- 
metric flow rate. 


Location of the controlling diffusional resistance 


The fraction of the total diffusional resistance 
which is located in the liquid lead phase can be 
expressed by rearrangement of equation (8). 

Fraction of Total Resistance in Liquid Lead 
Phase 

% 
(G,, Ly) Hi, Noa) [Pzn YZn Prly 4¥.)| dy 
—_ (*) 
Hog 


a 


As in the evaluation of the number of transfer 
units, an exponential decrease in the gas phase 
zine concentration from the bottom to the top 
of the column and a linear relationship between 
lead temperature and position in the column 
were assumed to evaluate the integral. 

The estimation of H, for zine absorption in 
molten lead required the application of the 
empirical correlation methods, discussed in the 
appendix, to the data derived from the carbon 
dioxide absorptions in water. The diffusivity of 
zine in molten the 
Srokes-EINsTEIN equation with the ionic radius 


lead was estimated using 
as the size of the diffusing species. The zine ionic 
radius was used because the published diffusion 
data of various metals in mercury support the 
concept of the movement of metals free of their 
conductance electrons in the liquid state. 
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For example in Run 2D: 
41, 

: P2n YZn - dy 

| Pry — %) 


Yo 
and (H,)z, — pp = 0-115 ft 


= 1-104 by graphical integration 


Fraction of total resistance in liquid lead 


phase = 0-0041. 


Hence, for Run 2D, only 0-4 per cent of the 
total resistance to mass transfer is located in 
the liquid lead phase. Similar calculations on the 
other runs indicated that the liquid lead phase 
were negligible, under the conditions 
of the zine absorptions carried out in this investi- 


resistances 


gation. 

It is important to note that the experimental 
conditions determine the contribution which 
the liquid lead phase makes to the over-all mass 
Thus, for example, if the 
experiments were carried out at very high gas 


transfer resistance. 


rates, very low liquid lead rates and much higher 
temperatures inspection of the equation (9) 
indicates that the fractional resistance of the 
lead phase could then become appreciable. 
For the experimental conditions of this investi- 
gation, the absorption of zinc vapour in molten 
lead is clearly controlled by the diffusional 
resistance in the gas phase, and any resistance 
in the liquid lead phase can be neglected. 
Under these conditions equation (8) reduces to 


Hog = He (10) 


COMPARISON OF ZINC AND AMMONIA 


ABSORPTION 


A comparison of the heights of gas film transfer 
units for zine vapour absorption in molten lead 
with those for ammonia absorption in water 


permitted an evaluation of the gas phase 
empirical correlation methods for predicting 
absorption rates of metallic vapours in liquid 
metals at elevated temperatures, an application 
far beyond the hitherto restricted range of such 
correlations. 

As a preliminary analysis of the zinc absorption 
data revealed anomalies when the inlet zine partial 
pressure was less than 1 mm Hg, such data are 
excluded from the correlation of zinc and ammonia 
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absorptions and the reasons for the omission 
discussed later. 

One packed height was used in the zinc absorp- 
tion experiments and the results must be con- 
sidered as apparent values of Hy, or Hg. It 
seems reasonable to assume that the end effects 
associated with the zinc absorptions are approxi- 
mately equal to those measured in the room 
temperature analogue. Therefore, these apparent 
values were converted to true heights of transfer 
units by assuming the end effects were equivalent 
to another 48 in. of packed height, before the 
comparison with the ammonia absorption data 
was made. 

In ammonia absorption, the dynamic hold-up 
is the effective component of the total liquid in 
the bed during the absorption. Therefore, as a 
logical start to the comparison of the two absorp- 
tion processes the conventional form of correla- 
tion given by equation (13) was used, with the 
dynamic hold-up being considered as the effective 
component in both cases. The exponent on the 
effective hold-up ratio was previously determined 
as 0-42. 

Inspection of the data plotted as HM, Se on 
Hy * vs. v, 9 « on a log-log graph showed 
that all the zine absorption runs up to and 
including Run 3B fell significantly the 
ammonia data, whilst the points for runs after 


below 


Run 3B were displaced above. The significance 
of this was not at first realized. It was only after 
a search of the experimental records, that the 
cause of the displacement of the data was believed 
identified. At the 
Run 3B, the stainless steel inlet analysis sampling 


to have been conclusion of 
line broke due to a combination of gross over- 
heating and zine attack with the result that the 
static hold-up of liquid lead in the column at the 
conclusion of Run 3B was exposed to oxidiz- 
conditions for about 15 and became 


ng min 


drossed up. It is not envisaged, however, that 
the liquid lead was completely converted to dross 
but merely drossed on the surface. Had the 
static the 


volume of hold-up would have changed con- 


hold-up been completely drossed, 
siderably, but hold-up estimations showed that 
the static hold-up volume was not significantly 


different after Run 3B than it was in the earlier 


runs once the packing had all its static hold-up 
sites occupied (i.e. after Run 4A). 

The significant feature, which determines the 
effectiveness of the semi-stagnant static hold-up 
in an absorption process, is the rate of movement 
of this component through the packing. For 
water irrigation of the packing, the rate of move- 
ment is very slow and hence equilibrium in the 
static hold-up is reached and it is then no longer 
However, if the semi- 
through the 
packing at a rate only slightly less than the 
recognised (Le. dynamic 
hold-up), then the static hold-up could be equally 


effective for absorption. 


stagnant = static moves 


hold-up 


flowing component 
as effective as the dynamic hold-up. 

It is suggested that the rate of movement of 
a liquid metal or similar dense medium static 
hold-up in a packed column would be considerably 
greater than that which has been observed for 
water. The higher kinetic energy of the flowing 
component could have the effect of “knocking-on” 
or displacing the semi-stagnant liquid through 
the sites of static hold-up, thereby increasing the 
over-all rate of movement of the so-called static 
hold-up under irrigated conditions. Evidence 
of this increased mobility was given previously 
when the observed effect of liquid rate on static 
hold-up of mercury was discussed by the author 
[6], whereas the static hold-up is unaffected by 
liquid rate when water is the irrigating phase. 

It is therefore tentatively proposed, pending 
that 
effective interfacial area for liquid metals irrigat- 
ing a packing be modified to take into account the 


further experimentation, the concept of 


increased mobility of the semi-stagnant liquid and 
the so-called static hold-up considered effective 
in an absorption process involving the flow of a 
liquid metal or similar medium. It is likewise 
considered that the of the 
liquid metal also facilitates an increased mobility 
If the free 
flow of the static hold-up was impaired by dross 
formation or similar effects, the mobility would be 
seriously affected and the movement would most 
likely be completely arrested or slowed down to 
such an extent that equilibrium would be reached 
and the surface rendered ineffective for further 
absorption, 


non-wetting flow 


of the semi-stagnant liquid hold-up. 
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A qualitative test was run to check the assump- 
tion that dross formation would prevent free 
the 
A quantity of steel rings was irrigated with molten 


movement of lead retained by the rings. 
lead and then partially oxidized. The rings were 
then quickly shaken before the lead froze. An 
examination of the rings, after such treatment, 
indicated that the drossed lead droplets were 
retained by the packing. A similar test carried 
out on clean mercury-irrigated packing showed 
that the mercury could be readily shaken from the 
packing. 

In view of the preceding discussion, it was 
therefore assumed that the total hold-up of liquid 
lead up to and including Run 3B was effective for 
After Run 3B only 
hold-up was considered to be effective for absorp- 


absorption. the dynamic 
tion, as the movement of the static hold-up would 
have been arrested by dross formation. 

The the 
making appropriate 


after 
the 


correlation of data obtained 


the 
effective hold-up ratio is shown in Fig. 8. 


amendment to 


The displacement of the zine absorption data 
above the ammonia data in Fig. 8 could possibly 
be due to a deficiency in the estimated values of 
Schmidt number or effective hold-up ratio, but 
it also suggests that the Reynolds number, based 
on the gas velocity, is not a good correlating 
modulus, as the ammonia data of these studies 
are particularly sensitive to a relative velocity 
correction, 

As a check to see if the gas inertia is the better 
correlating modulus, as suggested by equation (15), 
the data were replotted in Fig. 9 using the same 
ordinate as in Fig. 8 but with (v7 p)®* as the 
abscissa. The square root of the gas inertia was 
used in this plot so as to retain the same gas 
velocity dependency as in the previous graph. 

The vertical displacement of the two sets of 
data has been overcome using this method of 
correlation, and the data in Fig. 9 can be well 
represented by a single straight line. 

As the newly defined “relative jp” factor, 
developed by the author [7] from gas phase mass 
transfer studies with a disk column, is believed to 
have general application to all mass transfer 
studies, the new type of correlation proposed by 
equation (16), in which the relative interfacial 
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. 8. Correlation of zine and ammonia absorption data. 
@ Zine absorption in lead. 
O Ammonia absorption in water. 


velocity is considered the significant velocity 
term in the transfer process, was applied to the 
zinc and ammonia absorption data. The relative 
interfacial velocity used in the correlation of the 
data shown in Fig. 10 was taken as the sum of the 


gas and liquid phase average linear velocities. 
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Fic. 9. Correlation of zine and ammonia absorption 
@ Zine absorption in lead. 


© Amonnia absorption in water. 
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Fic. 10. Correlation of zine and ammonia absorption data. 


@ Zine absorption in lead. © Ammonia absorption in water. 


Although the experimental points in Fig. 10 — the conventional type of correlation using the gas 


exhibit a degree of scatter, the data can be velocity as the significant velocity term, as can 
represented by the single full line drawn on the be seen by comparing Fig. 8 with Fig. 10. 
graph. This is in marked contrast to the obvious As indicated previously, the zine absorption 


displacement of the two sets of data obtained for runs, in which the inlet zine partial pressure was 


06r- 











Iniet solute portico! pressure , mmig 


Fic. 11. The effect of inlet solute partial pressure. 


@ Zine absorption in lead, © Ammonia absorption in water. 
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The absorption of zine vapour in molten lead, 


less than Imm Hg, were excluded from the 
preceding analysis. The 


of these very low zinc inlet partial pressure runs 


anomalous behaviour 


can be seen by selecting the correlation form, 
observed as fitting the other data with the least 
seatter, and replotting the data in such a form 
as to show the effect of inlet solute partial pressure. 

Thus Fig. 11 is a plot of Hg (Se) °* (H,)** 
(uz py?’ partial pressure. 
Theory predicts that there would be no variation 
of the group plotted with the inlet partial pressure 


versus solute 


of the diffusing species. 

As can be seen from Fig. 11, there is clearly 
some extraneous factor coming into the analysis, 
when the zine inlet partial pressures are less than 
about I mm Hg. At these very low zine concen- 
trations, errors in the analytical procedures 
could be expected but these can hardly account 
for the very pronounced and consistent effect 
observed. The only feasible explanation appears 
to be that the low zine partial pressures were not 
effective in removing the oxygen and carbon 
dioxide from the gas circulation system before 
the runs were commenced, 

If there were no leaks in the system, any oxygen 
contamination would have been eliminated by 
reaction with zine vapour well before the runs 
were commenced, the was 
recycled, Thereafter, the only oxygen which 
could enter the system would be in the small 
commercial grade nitrogen make-up to compensate 
for leaks in the system, as the entire system was 


because nitrogen 


maintained above atmospheric pressure. As this 
nitrogen passed through a bed of granulated 
charcoal maintained at 1,000 C before coming into 
contact with the zinc vapour in the boiler, the 
effects of oxygen or carbon dioxide contamination 
could normally be considered as being eliminated. 

However, in the low zinc partial pressure runs, 


it is considered that the quantitative removal of 


oxygen and the conversion of carbon dioxide to 
carbon monoxide was not completed before the 
runs were commenced because of insuflicient zinc 
being present. Under these conditions, part of 
the zine entering the absorption column would be 
as zine oxide, which would pass through the 
column in particulate form without being 
absorbed by the molten lead. As the gas phase 


Countercurrent absorption in «a packed column 
analysis methods used in the _ investigation 
could not detect any difference between zinc oxide 
and zinc vapour, the gas phase analyses under 
these very low inlet zine partial pressure condi- 
tions would not be truly representative of the 
actual zine vapour present. 


CONCLUSION 


These studies clearly show the absorption of zine 
vapour in molten lead to be a diffusion controlled 
mass transfer operation, which is fundamentally 
the gas absorption 
carried out at room temperature with aqueous 


same as a_ conventional 
or organic liquids as the absorbing liquid phase. 
The contacting of zine vapour with molten 
lead, in such processes as the recently developed 
zine blast furnace, should be considered as an 
absorption not a condensation and the units used 
should be termed absorbers not condensers. The 
use of the nomenclature existing at present in 
metallurgical industry, in which the zine vapour 
molten valled a 
condenser, can lead to erroneous conclusions being 
drawn with respect to the potential performance 
of such units. These units may be operated 
with neither the lead or gas temperatures below 
the dew point of the zinc vapour as can be seen 


lead contacting chamber is 


by noting the temperature conditions of Runs 
2D and 6D, for example. In these runs, the lead 
temperatures were at no stage below the dew 
point of the inlet zinc vapour. The dew points 
for Runs 2D 6D 573°C and 576°C, 
respectively, whereas the lowest lead temperature 
in these runs were 610°C and 635°C, respectively. 
Such performance obviously could not be expected 


and are 


of a condenser. 

The most important effect of lead temperature 
is in determining the equilibrium zine partial 
pressure above the alloy formed by absorption. 
If this is less than the bulk gas phase zine partial 
pressure, diffusion of zinc from the gas phase to 
the liquid lead phase can be expected and, if 
suflicient time is allowed, the concentration in the 
gas phase will be reduced until the equilibrium 
condition for the particular lead surface composi- 
tion and temperature are reached, even though 
the gas is at no stage cooled to the dew point of 
its zinc vapour content, 
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Besides the equilibrium restriction in determin- 
ing the direction of the mass transfer, the other 
factor which limits the temperature of the lead 
phase is that increased lead temperatures result in 
increased vaporization of the liquid lead, although 
even operation at temperatures approaching 
1,000°C, the losses of lead would be very small. 
Imm Hg). 


Of the types of correlations discussed for gas 


(Vapour pressure of lead at 987°C 


phase mass transfer, the form based on the gas 
inertia or the relative interfacial velocity gave the 
best correlation of the zine absorption in molten 
lead and the ammonia absorption in water data. 
The correlations using the gas inertia appears 
slightly superior, but the relative merits of the 
two methods cannot be definitely established 
until further experimental work is undertaken. 
Despite the slightly better correlation obtained 
when the gas inertia is used, it is considered that 
the form of correlation using the relative interfacial 
velocity is fundamentally the best approach to the 
examination of a gas phase mass transfer in a 
packed column. 

The mobility of the semi-stagnant liquid hold- 
up was of considerable importance in assessing 
the effective hold-up of liquid lead used in the 
interpretation of the zinc absorption data. It was 
suggested that the total liquid lead hold-up was 
effective prior to the static hold-up becoming 
partially drossed, and thereafter only the dynamic 
hold-up was considered effective in the absorption 
process, 

Although the absorptions of zinc vapour in 
molten lead were controlled by the diffusional 
resistance of the gas phase under the conditions 
the 


liquid lead phase to the over-all mass transfer 


of this investigation, the contribution of 


resistance depends entirely on the conditions 
under which the absorption is carried out, as can 
be seen by inspection of the fundamental equation 
derived : 

¥) 

_ GM, 

Hy + — 
Ly, N OG « 
Mg 


om YZn - da 
I rly Ye) 


Hog 


If the value of the second term on the R.HLS. 
of the above equation had been significant the 
interpretation of the zinc absorption data would 


have been very much more diflicult, particularly 
as the values of H, should strictly also be included 
within the integral, because of the variation of 
liquid phase Schmidt number associated with the 
different inlet and outlet lead temperatures, 

Throughout this investigation the principal 
interest has been directed towards the analysis 
of gas phase controlled absorptions. For the 
absorption of a metallic vapour in a liquid metal 
at an elevated temperature, in which the principal 
diffusional resistance is located in the liquid phase, 
further research is definitely indicated. 

The application of room temperature analogues, 
for investigating experimentally difficult high 
temperature operations, has been demonstrated 
in these studies. The zine absorption in molten 
lead and the absorption of ammonia in water can 
be directly compared using appropriate empirical 
contain the 


correlations suitably modified to 


essential parameters of relative performance, 


APPENDIX 

Emrpinicat Mass TRANSFER CORRELATIONS 
(1) Correlation of liquid phase data 

Surrwoop and Hottoway [8] proposed that packed 

column liquid phase data be correlated according to 
equation (11): 

L " s 

HW, | “| 

ia 


PD 


(11) 


SHerwoop and Ho tvoway investigated a wide range 
of systems, with the result that the exponent (s) on the 
Schmidt number was found to be 0-53 in equation (11). 

The values of the exponent (#) and the constant («) 
depend on the size and type of packing, the nature of the 
surface of the packing and the type of flow, whether 
wetting or non-wetting, and therefore must be determined 
experimentally on the type and size of packing and flow 
for which the correlation is being developed. 


(2) Correlation of gas phase data 

Yosuipa [9| chose to correlate his data, for the vaporiza- 
tion of water by a series of gases in a packed column, by a 
form similar to that used by SHerwoop and Hottoway 
for liquid phase data. Equivalent forms in terms of the 
mass transfer coeflicients have been used by many investiga- 
tors of gas phase mass transfer. Moditication of these forms 
to include the true rather than superficial gas velocity 
gives equation (12), which may be considered as the 
conventional method of correlating gas phase data. 


b c 
Hg aus (%*) (*) 
4 PD 
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The use of the superficial liquid phase mass velocity, 
to account for the observed effect of liquid rate on gas 
phase mass transfer, is considered to be the principal 
weakness of these correlations, when examined from the 
aspect of application to the zine absorption studies of this 
investigation. If the correlation is developed by experi- 
mentation with aqueous systems, the extension of the 
correlation to a system, in which the liquid phase is 
molten lead, would be extremely questionable. Therefore, 
rather than using a liquid rate in the correlation, a term 
indicating the actual volume of liquid present in the pack- 
ing, which is effective for absorption, should be used in the 
correlation, It is proposed that the dimensionless quantity 
“the effective hold-up ratio” Hp be introduced for this 
purpose. The “ effective hold-up ratio” to be defined as 
the volume of liquid in the packing, which is effective in 
absorption, divided by the volume of the packed bed. 
By the introduction of this term, it is considered that a 
compromise has been reached between the use of effective 
interfacial area as suggested by SuuLMan ef al. (10, 11, 12) 
and the more readily achieved procedure of evaluating 
transfer rates per unit packed bed volume. Equation (12) 
can then be moditied to the new form given in equation (13). 


» b c 
Mg = 8 | e*) (“) To 
B PD 


Lyxcu and WinKt 


(15) 


[13] proposed that the rate of mass 
transfer may be a function of the inertia of the gas stream 
(v,? p) at high flow rates rather than the Reynolds number. 
The form of correlation suggested is given in equation 14. 


Hg = ¢ le p)s | ; | 


PD 


(14) 


Equation (14) was proposed tentatively pending further 
study. 
their data by such a method, their data can also be 
correlated conventionally, if the effects of interfacial 
velocity are taken into account, as demonstrated previously 
by the author [6]. 
of equation (14) can be made until further experimental 


Although these workers successfully correlated 


As no real assessment of the merits 


work is undertaken, a modified form involving the newly 
proposed effective hold-up ratio, as shown in equation 
(15), was also considered as a possible basis for the correla- 
tion of the zine absorption data. 


c 
Hg = B (vg? pf (*) HA (15) 
PD: 

The exponent of the Schmidt number in equations (12) 
to (15) is considered to be 4. The investigations supporting 
this exponent are very numerous. Although the recent 
studies of Lyncu and WiLke [12] and Yosuipa [8] obtained 
different exponents on the Schmidt number, these are 
believed to be due to interfacial velocity effects masking 
the true effect of the diffusivity, as was proposed previously 
by the author [6). 


Countercurrent absorption in a packed column 


The author's gas phase mass transfer studies [6] resulted 
in the introduction of a newly defined mass transfer 
factor, * the relative jp factor,” in which the significant 
velocity term in the mass transfer process is taken to be 
the relative interfacial velocity, rather than the true gas 
velocity. As the use of this new approach is considered 
to have general application, a possible form of correlating 
the zine absorption data of this investigation would be an 


equation derived from the * relative jp factor” as given 


by equation (16). 


: f 
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NOMENCLATURE 


diffusion coeflicient 
superfical molar mass velocity 
of the gas phase Ib mole /see ft* 
dynamic hold-up ratio 
effective hold-up ratio 
height of a gas film transfer unit ft 
height of an over-all gas phase transfer unit ft 
height of a liquid film transfer unit ft 
over-all volumetric gas phase coeflicient 
Ib mole ‘sec ft® atm 
volumetric gas film coeflicient 
Ib mole ‘sec ft® atm 
volumetric liquid film coefficient 
Ib mole sec ft® unit AC 
superficial mass velocity of the liquid phase 
Ib ‘see ft® 
superficial molar mass velocity of the liquid 
phase Ib mole ‘sec ft® 
number of over-all gas phase transfer units 
total pressure in the system atm 
= total pressure on system mm Hg 
equilibrium vapour pressure of pure zinc 
mm Hy 
partial pressure of zine in gas phase mm Hyg 
- Schmidt number = (y pD) 
= true gas velocity 


dimensionless 
ft see 
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relative velocity of the gas and liquid streams fluid density 
ft sec 

mole fraction of zine in liquid lead 

mole fraction of zine in the gas phase 


height of packed section of absorption column ft Subscripts 


. doe. fn. 8 = exponents in empirical equations 


. B. B’. B constants in empirical equations ; 
activity coefficient of zine in liquid lead-zine 2 = dilute end of countercurrent system 


alloy (standard state pure liquid zinc at the equilibrium value 


same temperature) .M. = logarithmic mean 
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Conserved-property diagrams for rate-process calculations—I 


Comparison of known diagrams 
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Abstract—The paper compares three known methods of representing the enthalpy and composition 
of binary mixtures, namely those in which these properties are referred respectively to the mass 
of mixture, the number of moles of mixture, and the mass of one component of the mixture. The 
water-steam-—air system is used as an example. Each of the diagrams can be used for calculating 
the rates of simultaneous transfer of heat and matter across a phase boundary. The graphical 
constructions necessary for determining the rates are described and their convenience discussed. 


The paper is primarily an assembly of known facts, prepared as a starting-point for the 
demonstration, to follow in Part Il, that the three charts are members of a much wider family, 
other members of which may be used with advantage in heat and mass transfer calculations. 


Résumé—Ce meémoire compare trois méthodes connues de représentation de lenthalpie et de 
la composition des mélanges binaires c’est-a-dire celles ol ces propriétés se rapportent respective- 
ment a la masse du mélange, au nombre de moles du mélange et a la masse de bun des composants 
du mélange. Le systéme eau-vapeur-air est utilisé comme exemple. Chacun des diagrammes 
peut étre utilisé pour le calcul des vitesses du transfert simultané de chaleur et de maticre a 
travers la limite des phases. Les constructions graphiques nécessaires a la determination des 
vitesses sont décrites et leurs avantages discuteés. 

Le mémoire compare d’abord un ensemble de faits connus servant de point de départ pour 
démontrer (& suivre dans la part 11), que les trois diagrammes font partie d'une famille trés vaste 
qui peut étre utilisée avec avantage dans les calculs de transfert de chaleur et de masse. 


Zusammenfassung —Diese Arbeit vergleicht drei bekannte Methoden zur Darstellung der 
Enthalpie und der Zusammensetzung von Zweistoffgemischen, insbesondere solche, in denen 
diese Groéssen entweder auf die Masse der Mischung, die Molzahl der Mischung oder die Masse 
einer Komponente der Mischung bezogen sind. Das System Wasser Dampf— Luft wird als Beispiel 
verwendet, Jedes der Diagramme kann zur Berechnung von gleichzeitigem Wirme — und 
Stofftransport durch eine Phasengranzfliche verwendet werden, Die graphischen Konstruktionen 
zur Bestimmung der UCbertragungsgeschwindigkeiten werden beschricben und ihre Vorteile 
diskutiert. 

Die Arbeit teilt im wesentlichen bekannte Tatsachen mit dient als Ausgangspunkt fiir den 
folvenden Teil Il. Hier wird gezeigt, dass drei Diagrammme Glieder einer viel grOsseren Schar von 
Diagrammen sind ; andere Glieder dieser Schar kOnnen mit Vorteil fiir die Berechnungen der 


Wiirme-und Stoffiibertragung verwendet werden. 


l INTRODUCTION binary mixtures. This has been used in one form 
or another by Poncuon [1], Savanrir [2], Monier 
[3], Keenan [4] and, perhaps most extensively, 
BoSxsakovie [5]. 


1.1 Representation of Thermodynamic Properties 
Tue thermodynamic properties of binary and 
ternary mixtures are often represented by means 
of charts, the ordinate and abscissa of which are (b) The triangular diagram for ternary 
properties featuring in conservation laws. Exam- mixtures. This is often used for delineating the 


ples are : phase boundaries, and representing the isotherms, 
(a) The enthalpy-composition diagram for of mixtures encountered in solvent-extraction 
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processes (see for example Surrwoop and Picrorp 
[6], or TreyBat [7)). 

(c) The ternary mixture 
re-plotted in the so-called “ solvent-free co- 
The then similar in 


data are often 


ordinates.” charts are 


appearance to enthalpy-composition charts of the 


PONCHON type. 

Recently, enthalpy-composition diagrams have 
been used, by BoSxsakovid [8, 9], Powerit [10], 
and Spatpine and Ty.er [11], to represent the 
properties of equilibrium mixtures of substances 
which undergo chemical reaction when the mixture 
ratio changes. These are particularly important 
in combustion technology. 


The mixing rule 

All these diagrams have an important common 
feature: they obey the “ Lever Rule ~ of mixing. 
This is illustrated in Fig. 1. It states that where 
two systems, each being a mixture of the two or 
three substances in question, are mixed adiabati- 
cally and at constant pressure, the state of the 
resulting mixture is represented by a point, Y, 
on the property diagram, which lies on the straight 
line joining the state-points X and Z of the initial 
systems. Further, Y divides XZ in such a ratio 
that 

XY — quantity of material in system Z 


(1) 


YZ quantity of material in system X 





Conserved property ——__——_> 


This means that Y lies at the correct point for 
the fulerum of a lever (Fig. 1). 

The “ Lever Rule” holds because conserved 
properties (enthalpy, fractional concentration) 
are used as co-ordinates ; it permits calculations 
involving material balances and the First Law of 
Thermodynamics to be represented, and indeed 
carried out, by simple graphical constructions. 

The term “quantity of material,’ in the 
“ Lever Rule,” has to be interpreted differently 
according to the actual co-ordinates used. Some- 


times it is measured in mass of mixture, sometimes 
in moles of mixture, and sometimes it refers to 
just a single one of the components. 

If ternary mixtures are in question, the restric- 
tion to adiabatic constant-pressure processes in 
the * Lever Rule ” can be dropped ; then however 
only the composition of the final mixture is 
determined, not its temperature. 


1.2 Representation of Rate Processes 


Enthalpy-composition charts have also been 
used for the separate task of calculating the rate 
at which heat and mass transfer will occur when 
mixtures of different phases, compositions and 
temperatures are in contact. This was first done, 
it appears, by BuseMann [12]; the method has 
been extended recently by BoSnsakovic, who has 
applied it to water-steam-air mixtures [13, 14], 








Conserved property ———> 
Fic, 1. 


Illustrating the “ Lever Rule.” 


184 





Conserved-property diagrams for rate-process calculations—lI. 


to chemical reactions catalysed on surfaces [8] 
and to the industrial gas-producer [9]. 

The method has the advantage, as compared 
with purely numerical methods, of introducing 
the non-linear inter-phase equilibrium relations, 
which provide boundary conditions in the rate- 
process calculations, without requiring trial-and- 
error. Its graphical nature enables calculations 
to be performed quickly and with little mental 
effort, and with an accuracy which is quite sufli- 
cient in view of the uncertainty in current know- 
ledge of transfer coeflicients, and even of inter- 
facial areas, in chemical engineering equipment. 

Similar methods have also been developed 
independently by the present author and applied 
to water-steam-air mixtures [15] and to the 
combustion of liquid and solid fuels [11]. 

1.3 Purpose and Content of the Present Paper 

The methods of BuseEMANN, BoSNJAKOVIC and 
the present author exhibit both similarities and 
differences. The differences lie chiefly in the 
choice of co-ordinate system for the thermo- 
dynamic charts and in the details of the graphical 
calculations ; the basic assumptions are very 
similar however. 

The present paper has arisen out of the effort 
to compare the validity, convenience and accuracy 
of the various methods, Its main result is to show 
that a still wider variety of method exists ; great 
freedom of choice of co-ordinate system is available 
and can be exploited with advantage. This will 
be discussed in Part IT. 

For the sake of concreteness, discussion is here 
restricted to the water-steam-air system, which 
is familiar from studies of psychrometry, drying, 
cooling-towers and the like. This means that one 
co-ordinate of the chart will represent composition 
and the other enthalpy. However, the theory to 
be presented is quite general ; the substances may 
react chemically with each other to produce any 
number of other compounds ; the concentration 
of a third material component may replace 
enthalpy as one co-ordinate; and indeed it is 
not even necessary that thermodynamic equili- 
brium throughout the whole 
physical system (though it must do so at certain 


should _ prevail 


points). 


2. Turee Types OF 
ENTHALPY-COMPOSITION 
DiaGRamM COMPARED 
2.1. Description of Charts 

(a) ‘* Mass-of-mixture ” basis 

Fig. 2(a) shows an enthalpy-composition chart 
for water-steam-—air mixtures at 1 atm pressure. 
The horizontal co-ordinate, m, is the fractional 
mass concentration of H,O (liquid and/or vapour) 
in the mixture; its units are lb,, H,O, lb,, mixture ; 
m == 0 for pure air, and m = 1 for pure water or 
steam. 

The vertical co-ordinate, h, is the enthalpy in 
B.t.u/lb,, mixture. The enthalpy is defined 
arbitrarily as having the value zero for pure air 
at 32°F and for pure water at 32°F. The vapour- 
phase mixture is supposed to obey Dalton’s Law 
(though this is not necessary for the major part 
of the subsequent argument). Solubility of air 
in water is neglected. 

The body of the chart is divided into two 
regions by the curve marked “ saturation line ” 
above lies the steam-plus-air (vapour) region ; 
below lies the steam-plus-air-plus-water (mixed- 


phase) region. The two regions are traversed by 
a few representative lines of constant temperature 
(isotherms). The region below 32°F is omitted as 


being of minor interest. 

Charts such as Fig. 2(a) have been used by 
the present author [15] ; in that reference however 
0 line was on the right of the diagram ; 
= 0 line has been placed on the left in 


the m 
here the m 
order to conform with the two other diagrams 
which are to be discussed. 
(b) “ Mole-of-mixture ” basis 

Fig. 2(b) shows an enthalpy-composition chart 
for H,O - air at 1 atm pressure which differs from 
Fig. 2(a) in having its enthalpy and composition 
based on the mole unit. Diagrams of this type 
have been presented by BoSxsaxkovic [14]. The 
abscissa, M, measures the composition as the 
number of lb moles of H,O in the mixture divided 
by the sum of the number of lb moles of H,O and 
the number of lb moles of air in the mixture ; 
the units of M are abbreviated to “ moles H,O 
moles mixture,” for short. One lb mole of H,0O is 
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500°F 
4 


200°F 


Fig. 2. 
(a) “ Mass-of-mixture “ basis (hy) 


- 
rhe 


of air is 20 lb... 
the 


mixture divided by sum of numbers of moles of 


ib, ; and one Ib mok 


‘molecular weight” of mixture 


HO and air) thus depends on composition (see 

section 2.2, equation 6), 

M — 1 for pure H,O. 
The 


enthalpy in 


M 0 for pure air and 


the 
zero 


vertical co-ordinate, H. is likewise 
B.t.u_ tb The 


state and other assumptions are the same as for 


mole mixture. 
the chart of Fig. 2(a). 
The general appearance of the chart is similar 
» that of Pig. 2(a); but the gas-phase isotherms 
ure rather more nearly parallel, and there are 


other distortions. 


(c) “ Mass-of-air ~ basis 

Fig. 2(c) shows an enthalpy-composition chart, 
for H,O 
duced by Mouser [3]. Its characteristic feature 
is that the co-ordinate properties are based on 


air at 1 atm pressure, of the type intro- 


unit mass of air. Thus the abscissa, 7, measures 


the composition as the number of lb,, of H,O in 


the mixture divided by the number of lb, of air 
mixture. So # = 0 for pure air; and 


i * for pure steam or water. Only the range 


in the 


* Moles-of-mixture ~ basis («) 


(mass of 


Vv 
(b) 


Enthalpy-composition charts for H,O- air at 1 atm pressure, 


* Mass-of-air “ basis. 


Os 4 
the range, no finite diagram of this type can 


0-3 is covered by Fig. 2(c); but whatever 


represent mixtures of pure H,O. 
The 


enthalpy of the mixture in B.t.a Ib, air. 


the 
The 
zero-enthalpy conditions exhibit a difference from 
0 for air at 32°F, as 
This is 


vertical co-ordinate, J. measures 


the previous cases: J 
before, but for steam (not water) at 32°F. 
why the 32°F gas-phase isotherm is now horizontal. 


N.B. 
used ¢ for enthalpy. 


Moniier, and all subsequent commentators, have 
i differs from J in having the same 
zero condition as # and 77. Then, in order to obtain a 
diagram of convenient shape, Monier uses “ slanting 
co-ordinates,” i.e. plots i vertically upwards from a base 
line which is not at right angles to the lines of constant 
vr; the slope of the line is chosen to make the gas-phase 
isotherm for 32°F horizontal. This practice is therefore 
essentially the same as the use of J, with reetangular 
¢ -ordinates, exhibited by Fig. 2(¢); the explicit use of / 
seems to be the more straightforward practice. 


Fig. 2(c) shows the division into two regions by 
the saturation line which has already appeared in 
Figs. 2(a) and 2(b). But the three charts are by 
no means geometrically similar to each other. 
What is their relationship ? 
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+) ed 
= 


Relations between the Co-ordinates 


We have now introduced three ways of des- 
cribing the state of a given H,O ~ air mixture : 
the h—m, the H—M, and the / 
systems. We shall later want to use the algebraic 


relations between these quantities which are 


x co-ordinate 


implied in their definitions. These are easily 


shown to be: 
For composition 
max (1 +.) 
' M 
(29 18) + [1 —(29,18)] M 
together with equivalent relations, such as 
r=m (1 —m) 


m 


aie (18, 29) + [1 


(18, 29) | m m= 


while the * molecular weight ” of the mixture, ,, 
is given by 

p = 18M m 

29 (1 M) (1 m) 


29 + (18 — 29) M 





For enthalpy 
h=-H pw 
i (1 +2) 


v hy, 
1+a7 


(9) 
where hy, = latent heat of vaporization of water 
at 32°F, 
such as 


together with equivalent expressions 


H h (10) 


l=i-wrvh, 


h—mh, ti (11) 
1 m 


The * Lever Rule” on the Three Charts 


The conservation-of-matter principle, and the 
Steady-Flow Energy Equation, can be written 


2.3 


for adiabatic constant-pressure mixing processes 
using any one of the three pairs of properties : 
y—-m, H-morl—x, Application of these principles 
to the mixing of systems X and Z to make a 
mixture Y, as discussed above, readily shows that 
the “ Lever Rule,” equation (1), is valid provided 
“ quantity “ is interpreted as mass of mixture in 
case (a) (h- m), as moles of mixture in case (b) 
(71 — M) and as mass of air in case (c) (1 — 2x). 

It can then be deduced that the position of the 
point Y on the line XZ of Fig. 1 is related to the 
ratio of the mass of system Z, w., to that of 


system Z, w,, by: 
(a) (h—m chart) 
XY Wz 
YZ Wx 
(b) (H—M chart) 
XY Wz 
YZ Wx 
(c) (1-2 chart) 
XY wy a 
YZ wz (1 


rN) 


14 
vz) \ ) 


where the subscripts X and Z indicate the system 
for 
evaluated. 


which the properties »~ and x should be 


With these interpretations, each of the charts 
can be used for the construction of the “ polar 
etc.), 
facilitate application of the conservation laws to 


diagrams (5, 4, 16, which so greatly 
chemical engineering equipment. These are too 
well known to be described here. 


2.4 Rate Process Calculations on the Three Charts 

In this section some of the more important 
graphical constructions will be described which 
heat transfer 
processes to be calculated. We consider an element 
of the exchange surface as illustrated in Fig. 3; 
this shows a water film in contact, through a 


enable simultaneous and mass 


thermal resistance, with a reservoir at temperature 
to; the free surface of the liquid is swept by a 
steam-air mixture. This might represent con- 
ditions in a cooling-tower, in which case the 
reservoir is the bulk of the water stream at the 
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Gas stream 





Klement of exchange surface, with heat and 


mass fluxes. 


section in question, and the thermal resistance 
referred to above is the boundary layer of the 
water stream. A similar situation arises in con- 
densers, in which case the reservoir is the cooling 
fluid, while the thermal resistance is made up of 
the cooling-fluid boundary layer, the tube wall 
and its adher.ng scale and the water film on the 
outside. 

For the sake of simplicity, it will be assumed 
that radiative heat transfer between the water 
surface and its surroundings is negligible, and that 
the enthalpy increase of the water flowing in 
the film can be disregarded. The more general 
case is dealt with in [15]. 

We suppose that the state of the gas stream 
and the temperature of the reservoir are given ; 
in addition the gas-side heat transfer coectlicient 
and the over-all heat transfer coeflicient between 
the water film and the reservoir are supposed 
known. The problem is to determine, first of all, 
the temperature of the water surface and so the 
state of the gas immediately in contact with it ; 
thereafter the rate of mass transfer (condensation 
or vaporization), the rates of heat transfer from 
gas to water and from water to reservoir, and the 
direction of change of the gas state, are all to be 
determined. 

It is assumed that the liquid (L) and the gas 
(S), in contact at the phase boundary, are in 
thermodynamic equilibrium. 


(a) On the h—m chart 
The method developed by the present author 
[15, 11], utilizes the * mass-of-mixture ” basis for 


enthalpy and composition, The following con- 
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structions and formulae are extracted from [15] ; 
Their derivations are outlined in Appendix A. 
They are also valid for mixtures which undergo 
chemical reactions, 


Determination of the interface states 1. and S. 
(See Fig. 4). 
(i) Locate O on the line m — 1 corresponding 
to water at temperature {,). 
(ii) Locate H on the line m 1 such that 
the vertical distance OH (measured upwards, and 
in enthalpy units) is given by 


ey (te to) 


OH = 
a - [oo —1)](U ge) 


(15) 


(iii) Draw the line HG and tind A on it such 


that 


AG 
HA 


(16) 
(iv) Locate L and 8S, which lie at the extremi- 
ties of the mixed-phase isotherm through A. L 
lies on the line m 1: S lies on the saturation 
line. 
N.B. 
also, the procedure simplifies ; 
© so that (15) does not have to be evaluated: and the 


When the Lewis Number equals unity, and soo = 1 
for then H. coincides with 


second term of (16) disappears, 


Transfer rates (Pig. 4). 
G to intersect the gas-phase isotherm through S 
at D. Produce LS to cut DG at C. Plot K on DG 
such that 


Drop a vertical line from 


DK 
DG 


N.B. K and G coincide for ¢ 1. 


Then the mass transfer rate is given by 
"a CS 
(18) 
SL 


the heat transfer rate to the reservoir is given by 


qd’; o 


og 
= 


CK (19) 


and the heat transfer rate from the gas to the 
liquid surface is given by 


DG (20) 
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Change of gas condition (Fig. 5). Produce KS to 
cut the line m — Lat J. Let T be the intersection 
of GJ with the vertical through S. 

Then the state of the gas stream changes, as a 
result of contact with the liquid, in such a way 
that the corresponding state-point G moves 
towards J. 

More specifically, if the mass flow rate of the 
local gas stream is ri,, the elementary contact 
area is 5a, the state of the gas before contact is 
represented by G, and after contact by G,, then 
G, lies on the line G,J at a point such that 

GG, TJ gdba 
GJ TG, 


om 


PY a say (21) 


The expression on the right-hand side of (21) is 


the increment in the number of gas-side mass 
transfer units. 


(b) On the H—M chart 
[14] 

procedures permitting the calculation of heat and 
the *‘ mole-of-mixture ” 
basis ; his procedures differ from, and are some- 


BoSNJAKOVIC has described graphical 


mass transfer using 
what more complicated than, the constructions 
which will now be recommended for use with this 
chart ; but they are based on the same principles 


and assumptions, 


Satn line 


m\M) 


Fig, 4. 
(Erratum: DS should be prolonged to intersect OH at V) 


Rate-process constructions. 


Satn. line 


Fic. 5. Change of gas state. 

The constructions now to be described do not 
differ qualitatively from those of Figs. 4 and 5; 
the same diagrams can therefore be referred to. 
Proofs will not be given as yet ; they involve use 
of a result to be derived later (Part II, Section 
1.3). 


Determination of the interface state, Land S$ (Fig. 4). 

(i) As for case (a). 

(ii) As for case (a), except that the expression 
on the right-hand side of (15) must be multiplied 
by 18, the molecular weight of H,O. 

(iii) As for case (a), except that the expression 
on the right-hand side of (16) must be multiplied 
by w/18 where pe is the “* molecular weight ” 
corresponding to state G. It is convenient to 
provide a horizontal scale for this quantity on 
the diagram. 

(iv) As for case (a). 

Transfer rates (Fig. 4). D, C and K are points 
determined as in case (a). Then the transfer 
rates are calculated from the following equations, 
which parallel (18), (19) and (20) ; 


(22) 
(28) 


(24) 
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Change of gas condition (Fig. 5). The construe- 
tion and its evaluation are identical with those 
of case (a): equation (21) still holds. 


v chart 


(c) Onthel 
The early rate-process constructions developed 
by BuseMann [12] were made on the J — x chart. 
Much of the of BoSxsakovic [8, 13], 
including that on the carbon- air system [9], has 
The 
constructions now to be presented differ somewhat 
from his in the interests of simplicity and of 
conformity with those of cases (a) and (b). 
Determination of the interface states, L. and 5S 
(Fig. 6). The characteristic feature of the J —xr 
chart is that the pure-H,O line is at infinity. 


work 


also been done with charts of this type. 


Fic. 6. Rate-process and change of gas state. 
O, H, Land J can only be established as directions, 
not as points. The procedure therefore is : 


(i) Note the isotherm, tf, 
(ii) Evaluate ¢,,, from 


tg — lo 


‘n 1 + [oe (eo —1)](U ge) 


and note the corresponding isotherm. 
Two procedures are now possible. We first 
give one which is similar to those of cases (a) and 


(b): 


(iii) Draw a line through G parallel to the 
mixed-phase branch of the ¢,; isotherm. Find the 
point A on this line such that 


U 


£ ey 


o 4 (26) 
(iv) As for case (a). 

This procedure will fail if A lies beyond the 
right-hand edge of the available chart, as is 
likely tooceur when Ug cyislarge. The following 
alternative procedure may be used : 

(iii) Drop a vertical through G to intersect 
the gas-phase isotherm ¢,, (produced if necessary) 
at E. 
N.B. ME 


, Jf ° U | , 
MG / \1+——. . where M is the 
/ | o—1 ge} 
intersection of the vertical with the gas-phase isotherm /), 
produced if necessary. 


Find S such that 
CG aU 
ED gc 


(iv) 
+ (¢ —1) 


where C and D lie on the vertical through G and 
respectively on the mixed-phase and gas-phase 
isotherms through S. L, of course, is a point at 
infinity on the mixed-phase isotherm through S. 

The determination of S to satisfy (27) proceeds 
very quickly by trial-and-error: a point C is 
guessed ; D is immediately located and the ratio 
CG ED evaluated ; comparison with the pres- 
cribed value of the ratio shows immediately 
whether C should be moved up or down, 

This procedure can also be used, without 
modification, on the h —m or H — M charts. 

Transfer rates (Fig. 6). D and C have already 
been determined; K can be determined as for 
case (a). Then the transfer rates are calculated 
from the following equations, which parallel (18), 
(19) and (20), and (22), (23) and (24): 


CS h 
(28) 
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It will be noted that (1 
taken by « in the H — M formulation, 


x) occupies the place 


Direction of change of gas state. G, is related to 
G, as for case (a) except that, since the point J 
is at infinity, G,G, is drawn parallel to KS, For 
the same reason, equation (21), which is still 
valid, simplities to 


GG, 
GJ 


g-sa 
(31) 
og 
Comparison of the Merits of the Three Charts 
for Rate- Process Calculations 
The h 
The 


mally, and in the author's opinion conceptually 


m chart 


“mass-of-mixture”’ basis provides for- 
also, the simplest of the three sets of procedures. 
No y's or 1 + 2's appear in the equations used 
for evaluation. Only in respect of the evaluation 
of G,G, (via equation 21), does the /— m formula- 
tion involve a more complex expression than the 
I — x formulation. 

The reason for the simplicity of equations using 
the “ mass-of-mixture ” basis lies essentially in 
the fact that the tluid-mechanical processes which 
occur in the boundary layer adjacent the surface 
are responsive to the mass velocity of the stream 
rather than to the * mole velocity ~ or the ** mass- 
of-air” velocity. g can be evaluated without 
direct reference to molecular weight or compo- 
sition. 

The main disadvantage of the h—m chart is 
that its scale is often such that intricate construc- 
tions have to be performed on a small section of 
the chart. For H,O 


covering states near that of the atmosphere is 


air the region of the chart 


often too small for adequate accuracy of construc- 
tion. 


(b) The H— M chart 


The “ mole-of-mixture ” basis introduces the 
position-dependent quantity jz¢ into a number of 
formulae connected with the graphical construec- 
tions and their interpretation. This makes the 
significance of lines and ratios on the chart 
somewhat harder to grasp. 


There is a slight enlargement, as compared 


with the h— m chart, of the important atmospheric 
region for H,O But this is not usually 
suflicient, in the author's opinion, to compensate 


air. 


for the greater conceptual difficulty. The use of 
H—-M 
advisable. 

However it will later be shown that the important 
feature of the mole basis is that different mass 


co-ordinates does not therefore seem 


units are used for the two components. We shall 
see that there can be great advantage in using 
two different mass units, but that there is no 
reason why these units should be related tothe 


molecular weights of the substances. 


r chart 


(c) Thel 

In addition to having the position-dependent 
quantity (1 + 2) in many of the formulae, the 
meanings of the constructions on this chart are 
rendered diflicult to grasp by the fact that all 
pure-H,O points lie at infinity. This absence 
of the m = 1 line sometimes requires the use of 
slightly more involved constructions for determin- 
ing the interface state (either BoSnJakovic's, or 
those presented above) than are needed for the 
hm method. This is true not only of the rate- 


process constructions which are chiefly in question 


here, but also of the polar diagrams representing 
mass and energy balances. 

1 line absent, but the J— 
chart is incapable of representing mixtures which 


Not only is the m 


have an appreciable proportion of H,O in them. 
Heat and mass transfer in condensers, to name 
an important practical example, cannot be studied 
on the I 

However the J 
that it can easily be drawn to 


r chart. 

rv chart has the great advantage 
a suitable scale 
for the solution of problems involving relatively 
small moisture contents, such as occur in heating- 
and-ventilating practice and in cooling-towers. 
A minor advantage is that the formula for the 
change of gas condition, equation (31), is simpler 
than the corresponding one for the h- m and 
Ht — M charts. 


2.6 Summarizing Remarks 


It appears that there are striking similarities 
between the three charts discussed. Each has 


its own advantages and disadvantages, and which 
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is the best to use depends on the problem to be 
solved, whether it involves high or low moisture 
contents for example. 

Before trying to decide which chart is best used 
in various circumstances, we should inquire 
whether these three bases for the mixture quantity 
are the only ones which can be used. This is 
the subject of Part II. 


APPENDIX A, 
+ Derivation of Rate- Process Constructions 


A.1 Determination of interface states L. and S 


From the theory of heat and mass transfer (17), we have 


h Aca 
g. Ae = (A.1) 
(V's m”) Nyy 


and 
iW nm me 
.. (A.2) 
o ms 1 
Consideration of similar triangles in Fig. 4 shows that 
(A.2) can take the form : 
g be 
o LV 


m (A.3) 


From Fig. 4 and the fact that hyg LV, (A.1) becomes : 


DG 


g. (A.4 
* (q'p mh) + LN 


where V is the intersection of DS produced and the line 
m 1. 


From (A.4) and the definition of U', viz. 
Vs U (t, — to) 
(Uc cf) OL 
we deduce : 
q°DG (U 7) OL 
LV 


m 


Eliminating m” from (A.3) and (A.6), we find : 


CG = (Ue /¢g) OL — (¢ — 1) DG (A.7) 


Introducing into (A.7) the definition of point H (equation 
15) in the form : 
(cy ©) DG + OL 
1)] .(U ge) 


OH (A.8) 


+ [oe (oe 
we obtain : 
CG aU 
HL ge 


But, from similar triangles in Fig. 4, 


CG AG 
HL HA 


(A.10) 


from which the required result (16) follows immediately. 


Transfer rates 

Equation (18) follows from (A.2) by similar triangles. 
Equation (19) follows from (A.4), (18), the definition (17) 
and similar triangles. 
Equation (20) may be regarded as a definition of q",. 


Change of gas condition 


From (18), (19) and similar triangles, 


LJ V7 m” (A.11) 


Therefore the actual heat and mass fluxes through the 
surface are equivalent to a mass flux without heat flux 
It follows from the 
Lever Rule that the gas state moves towards J and that 


with the material entering at state J. 


GiG, om". ba 


(A.13) 
Ged mo 


Introduction of (A.2) and similar triangles thereupon 
lead directly to the required result (21). 


NOTATIONS 
State-point of gas stream 
State-point of liquid surface 
State-point of gas in contact with liquid surface 
Composition of G and temperature of S 
State-point of water at reservoir temperature fo 
Over-all 
surface to reservoir 


heat transfer coeflicient from liquid 


Gas-side heat transfer coeflicient divided by c 
(g Stanton 
gas stream ; 


Number times mass velocity in 
its units are Ib, ft® hr) 
Specific heat of mixture G at constant pressure 
Specific heat of water at constant pressure 


Mass transfer rate from liquid to gas phase 
Ib, ft? hr 


Heat transfer rate to reservoir (across lower 
surface of dotted control volume in Fig. 3) 
(across upper surface of dotted control volume 
in Fig. 3) 

Function of Lewis Number (Le = diffusion 
coeflicient of steam in air divided by thermal 
diffusivity of gas mixture) appearing in Chilton 
Colburn analogy between heat and mass transfer, 
where it is assumed o 


(1/2) and (2/3) [15]. 


(Ley-” when n is between 
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The longitudinal diffusion coefficient in flow through a tube with stagnant 
pockets 
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Abstract—In Turner's * Model | flow structure a fluid is assumed to flow through a tube whose 


wall leads to pockets of varying lengths. It is here shown that the presence of these pockets 
influences the Taylor diffusion coefficient for flow in the tube and a general method for calculating 
this is given 

Résume— Dans le systéme dccoulement présenté par Turner dans son “ Modéle 1° un thaick 
est suppose s‘ccouler a travers un tube dont la parol aboutit a des poches de longucurs variables. 
Lo auteur montre ici que la présence de ces poches influe sur le coeflicient de diffusion de Taylor 


pour Téecoulement dans le tube et il donne une méthode générale pour calculer ce coeflicient. 


Zusammenfassung 
durch cine Flissigkeit zustande, 
Linge besitzen. 
coeflizienten fur dic 
allvemeine Methode angegeben 


TURNER {1} model of tlow 


systems in which there are stagnant pockets of 


has pre yposed a 


various depths opening on the flow channels. 
The diffusion coefficient in these pockets is the 
molecular diffusion D whilst in the flowing channel 
a longitudinal diffusion coeflicient D is assumed. 
Neither in Turner's work nor in a later paper 
[2] on the same model is a value ascribed to D 
and one might be inclined to write for it the 
apparent longitudinal diffusion —_coeflicient 
D «tl? a@ D found by TAYLOR [3] (U 
radius of circle having the same 


mean 
velocity, a 
area, « is a shape factor, « = 1 48 for circular 
cross-section). This however is not correct and 
it is our purpose of this note to show how « 
depends not only on the cross-section and flow 
profile, as in Taytor’s analysis, but also on 8. 
The result is related to the behaviour of the tubular 
chromatographic column, which has been treated 
more generally elsewhere [4, 5], but cannot be 
immediately deduced from it. The result will be 
deduced from a comparison of two methods used 


to find the apparent diffusion coeflicient, the one 


Das Verhalten ciner StrOmung nach dem * 


kis wird hier vezeiet, dass die 
Rohrstromung nach Taylor beeinflusst ; 


Modell 1 von Turner kommt 


die durch ein Rohr stromt, dessen Winde Taschen verschiedener 


\nwesenheit dieser Taschen den Diffusicns- 


zu seiner Berechnung wird cine 


for Tayior’s analysis [6] and the other for 


Turner's model [7]. 


1 Two Meruonps or DETERMINING 


Apparent Dirruston CorrricieENtTs 


To see this comparison unencumbered consider 
a tube of length L through which a fluid is flowing 
with uniform speed U and in which a solute is 
If « be 


the concentration of solute it satisfies the equation 


diffusing with diffusion coeflicient D*. 


(1) 
where « is distance along the tube and ¢ the time. 
In the first method it is subject to the conditions 

st 


le U f(t), 4 


and we take the inlet concentration to be a unit 
amount of solute introduced at zero time, Le. 
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Cy (t) = 8(t). The emergent concentration as a 
function of time c (L, t) will be a peak representing 
the probability of a molecule having residence 
time ¢ and can be used to define a mean / and 
variance o*, where 


i [rou t) dt 


0 


o? | (t — i)? ¢(L, t) dt 


Van DER Laan [7] has shown how these quantities 


may be calculated from the Laplace transform of 


the solution of (1) and (2) without inversion. In 
this case 


2 D* L 


m+ O(L") 


whence 


b* 


Using this method on TUrNeR’s model it is shown 
elsewhere that, with a diffusion coeflicient D in 
the channel and a distribution of pockets such 
that 8(/)dl is the volume of pocket of length 
between / and | + dl per unit volume of channel, 

the apparent diffusion coeflicient is 

D OP fase 
= - ‘ 3 
3D (1 + py. p(t) (3) 
0 


There is however another way of looking at 
(1). very tube 

2 <«# < © and suppose that at time t = 0a 
unit amount of solute is concentrated on the plane 
a = 0. We may ask how the solute distributes 
itself as time goes on and regard the solution 
c(#,t) as a probability distribution giving the 


equation Consider a long 


probability density of a molecule reaching a 
point 2 in time ¢, From this viewpoint we would 
study the spatial moments 


«a 
. 


m, (t) = | re(a,t)dx 


* 
= @ 


a 


ms, (t) = 


7 
x 


(a — m,)* ¢ (a, t) da 


To do this it is convenient to take a moving 
origin with «=0 at t=0 and 
travelling with speed V, will be later 
chosen to be the speed of centre of mass of the 
— Vt then 


coincident 
which 


solute. Let z =a 


(4) 


Now let 


cl” (t) 2’ c(z, t) dz 


so c'”’ (0) 1, ce?’ (0) Then 


multiplying equation (4) by 2” and integrating 


0. Pp 


gives 


de' p) 


V) Apt) 
dt sit 


b* p(p 


As we should expect de‘’’ dt — 0 showing that 


the total amount of solute is constant. 
de) 
at 
and if V U, A” Thus the 
centre of mass travels with the uniform velocity 
of the stream. Finally putting p = 2 
de? 


dt 


is constant 0. 


2D* or c* 2D*t 


Thus 


m, = Vt= Ut, m, = 2D*t 


and 
D* = Um, 2m. (5) 


This technique has been applied to Tay or’s 
problem of simultaneous diffusion and convection 
by a non-uniform velocity, and will now be 
applied to TurNeER’s model to obtain a formula 
which on comparison with (3) gives D. 


2, APPLICATION TO TURNER’S MODEL 
We will consider circular channels of radius a 
but the method can be modified for any other 


cross-section without difliculty. On the wall of 
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this tube there open the mouths of the stagnant 
pores or pockets and it is assumed that the 
distribution of these does not vary along the 
length of the tube. 
Let 

8 (1) dl = volume of pockets of lengths 
between / and / + dl per unit 


volume of channel 


B = | B(l)dl = total volume of stagnant 


o 
channel 


¢, (@, y, O concentration of solute at 
depth y in a pocket of length 
l at the 


channel and time ¢ 


position r along 


concentration of solute in 
flowing channel at radius +, 
distance x and time f 

velocity of stream at radius 


r, U 


Ud(r) 


mean velocity. 


Then since the total area of pockets of lengths 
between 1 and / + dl is «(/) B(l) 1 per unit 
volume of channel we can, by attributing the wall 
area which is not a pocket mouth to pockets of 
length zero, ensure that 


x (I) dl | 8 (Il dl 
o o 
without changing the value of 8. The equations 
satisfied by ¢, and ¢ are thus 


\2, 
p° (6) 
dy? 


c (a, 2, (7) 


€, (a, 0, t) 


” a F dc 
=— fa(l)(—'}dl, +r 
a] (55 0 


u 


pocket per unit volume of 


If the moving co-ordinate system is introduced 
by zs = «2 — Jt equations (6) and (9) become 


(12) 


\ . \2n \c 
‘ . , », © 
dt )z? > 


(13) 


Now define the moments 


x 
. 


| 3? ¢ (z, r,t) dz = Cc?’ (r, ty 


. 
a 


| ~P C, (Y, Zz. t) dz ef?’ (y, t) 


pe”? (0) | 2rc (r,t) dr 


£ 
. . 


| x (1) | cP?’ (y, t) dy 


. . 
v0 “o 


Then multiplying through equations (12) and (13) 
and conditions (7), (8), (10), (11) by 2” 
integrating we have 


Dp (p 
; p(Ud 


(18) 
(19) 
Finally integrating with respect to r and y and 


using the boundary conditions (15), (16), (18) 
and (19) 
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Dp(p— 1) . | arc? *) dr 
a* 


(p) 


du 

dt 
pu at (p-l) y wie—t) » 

ob at | ’Urdirjec dr —pVp (20) 


0 


Initially the solute is concentrated in the plane 
z=—0 so ec” rd (1+ ,)y", pO ] 
GP = of”? = py‘?! -0O, att — 0. Putting 


0 so x‘? is constant, 


and 
0, p 

p = Oin (20) gives du’ dt 
which is to be expected since it is the total 
amwunt of solute present, Also equations (14)-(19) 


(0) ’ 


are satisfied if c°’ and ¢;"’ are constant. Thus for 


all t 
(21) 


For p = 1 equation (20) gives 


r a 
du! y ' 

a Dv» . adr 
dt (1 + Ba. 2rd (r) dr 


mean value of d& is 1. Thus if 


| B) the centre of mass is stationary 


since the 
V= 0/1 
in the moving co-ordinate system, being in fact 


at the moving origin. 


) 1) 


In seeking solutions for c'’ and ¢;'? it has been 
shown [6] that we need only concern ourselves 
with the part of the solution independent of time. 
The time dependent will ultimately grow to a 
small constant in the expression for «'*’, but since 
nw”? is increasing linearly with time such a constant 
is ultimately negligible. For the time independent 
part of the solution c;'’ equations (17)-(19) give 


de; 1) 


Vy 
ef) =” (a) paap 
id l 


Jaw [ «i dy dl = Be'"’ (a) 
0 0 
a ay l)P dl 
~ BD(1 4 5 | saci 


0 


Also the time independent part of c'’ satisfies 


D U 
r or\" ) + Py laa + B 


l 
P(r) - | 2r d(r) dr 
a 


so that 
P(v) = 0, D(a) 


and 
dct) laa 


yr ~ 2D(1 + 8) 


l d U 
b(r) 4 r 
r 2D(1 + BP 
Since @(r) 0 (r*) at r o, dc" Yr vanishes at 
r 0 and the condition (23) is identically satisfied 
atr — a. Integrating again 
U a 


Al) | 
oo) * spa 


c'(r) 


U , : 
wD aye" a*) 
and c''’ (a) can be determined by the condition 
(1) 

fe 0 


x 


| B(l) Fadl 


ce" (a) 


3D (1 B)* | 


a 
° 


_— » 1 r@Pi(r)dr (25) 
+ B)P ay | 

0 

Let p = 2 in equation (20) and substituting from 
(22) and (24) 

2« Ua 

D(1 + 8) 


a See ee 
+ spa cap} 2M! d ‘26) 
0 


dp 2D 
dt 14+ 8 


where 
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and p =r a. Thus 

m, = Ut/(1 + 8) 

Mm, = 2Kt 
where 2K is an abbreviation for the right-hand 
side of (26). Then by (5) 


Um, l?a* 


= Sm, =D+« i 


U? , 
+ spa 4 | B(l)P dl (28) 
0 


Comparing (28) with (3) we see that the correct 
value for D, the effective longitudinal diffusion 
coellicient in TurNeR’s model, 


D= D+ «U*%a D (29) 


where « is a constant depending not only on ¢ (r) 
but also on 8. 
For the parabolic protile ¢(r) = 2 (1 


1+ 68 +118? 
(1 + By 


This reduces to TayLor’s value when 8 = 0 but 
rises to eleven times this value when § is large. 
This is physically reasonable for the longitudinal 
dispersion is caused by the solute spending time in 


r? a*) 


slowly moving regions. Thus if there is a large 
region of zero speed, i.e. 8 is large, the longitudinal 
diffusion is enhanced. 


NOTATION 
radius of circular channel 
concentration in channel 


concentration in pocket of length / 
@ 


( 2P « dz 
ry 
molecular diffusion coefticient 


p™ moment of ¢ 


apparent diffusion coefficient in a flow system 
longitudinal diffusion coefficient in Turner's 
model 
length of pocket 
length of channel 
i*® spatial moment of concentration distribution 
time 
mean residence time 
mean flow velocity 
velocity of centre of mass of solute 
length along channel 
depth into pocket 
a Vit 
Bilt 
probability density of pockets of length / 
Dirac’s delta function 
shape factor 
p™ moment of total solute 
variance of residence time 
ratio of velocity at r to mean velocity 

; 

2rdéirjdr 


a. 
0 
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AMUNDSON 


( Received 7 April 1959) 
Abstract—The previous analysis of reactor stability and control is here extended to show the 
new features which arise when a derivative mode is added to the proportional controller. It is 
found that more care is needed to ensure full control of the whole phase plane, for the possibility 
of partial control is greater than in the previous example. It is also necessary to consider the 
effect of delay in the controller, These phenomena emphasize the care that must be taken in a 


satisfactory stability analysis. 


Résumé 


montrer les nouveaux caractéres qui se préesentent quand une action dérivée est ajoutée a Vaction 


Létude precedente sur la stabilite et le contréle d'un reacteur est poursuivie pour 
proportionnelle du contréleur. Les auteurs ont trouve que pour assurer un controle total du 
palier de stabilitée, il faut apporter plus de soin, la possibilité dun controle partiel étant plus 
loeffet de retard du 
Le phénomene souligne les précautions a prendre pour une analyse de stabilité satisfaisante. 


grande que dans Texemple precedent. controleur est aussi a considérer. 


Zusammenfassung —Kine friihere Untersuchung tber die Stabilitéit und Regelung eines 
Reaktors ist hier fiir den Fall erweitert, dass dem Proportionalregler ein Vorhalt gegeben wird. 
Es zeigt sich, dass man in diesem Falle stirker auf die volle Regelung im ganzen Phasenbereich 
achten muss, weil die Moglichkeit einer Teilregelung grosser ist als im friiheren Beispiel. Auch 
die VerzOgerung in dem Regler muss in Betracht gezogen werden, Auf alle diese Erscheinungen 


ist bei ciner ausreichenden Stebilitatsuntersuchunyg zu achten. 


THE PRECEDING papers of this series [1] have and a slightly different reduction to dimensionless 


outlined the methods of analysing the stability variables adopted, 
and control of a single stirred tank reactor. The 


possibility of control by a mixture of proportional 1. Tae SvsTex 


and derivative control was established (Part I, = ; ; 
anet ' Ihe system to be considered is a single stirred 
Section 4), but the evolution of the system only ’ 
This 


paper will show the new features which arise 


tank reactor of volume V in which a second-order 
If A is the 
concentration of the reactant which is being used 


fully described for proportional control. ‘ tae : 
: irreversible reaction is taking place. 


when a measure of derivative control is added to . : 
at a rate /A® per unit volume a mass balance 


the proportional controller, As might be expected 
from the very responsive nature of derivative 
control the possibility of 
It is also found necessary to consider 


hard oscillations is 
greater. 
the effect of a small delay in the controller, for 
however small this may be it imposes on the 
derivative control gain constant a restriction that 


A 


second-order reaction is used in this discussion 


could not otherwise have been anticipated. 


vives 
dA 


I = Q(Ay 


A 
de 


~kV AP (1.1) 
where @ is the time, Q the flow rate and A, the 
inlet concentration. / is the reaction rate constant, 
a function of temperature k = kh, exp(— E/RT). 

If the temperature is controlled by a cooling 
coil through which a coolant flows at a rate Q,, 
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entering at temperature 7 and leaving at 


T. + AT., a heat balance gives the equation 


AH) iAP 


(. P, c, AT. 


c, are specific heats of reactants and 


Qep (T, rT) 
(1.2) 
where ¢, 
coolant, and 
. p, are densities of reactants and coolant 
respectively, 
T, is the inlet reactant temperature, 
AH is the heat evolved for each mole of A 


destroyed. 


Since the mean temperature of the coolant ts 
T. + (1 2) AT, a heat balance on the cooling 
coil gives 


Q. Pe Ce AT, 
where A, is the area of cooling coil and / the heat 
By solving (1.3) for AT. 
and substituting into the last term of (1.2) we 


h A,(T — T (1 2) AT.) (1.3) 


c 


transfer coetlicient. 


have 
- r a 
Q. p. ¢. AT, pea T.) 
(Tt 
where 


F 2 Q. P. € h A, 


If h were constant F would be proportional to 
Q.. This would be the case if the resistance to 
heat transfer were wholly on the reactor side, for 


the value of h would then be set by the level of 


agitation. If the resistance were wholly on the 


inside of the cooling tube we should expect / 


with q., 


to vary a popular correlation gives 


Table 1. 


W, Tierney, R. Aris and N. R, Amunpson 


hh, Q**. The combined effect of these will 
give 


l I l 
ho hQ*? h, 


2 pete {Q"" | Q) 
A, | h, hy! 


i rd] (1.5) 


Thus F is a function of Q. and gives a dimension- 
less form of the flow rate, and eliminating Q 
between these two equations U can be expressed 


(1.5) 


whence 


then 


U'(q.) 2. Pe ©, [ 


”» . 1 ’ 
1 = Pe c Q. 
A, h, 


as a function of h. The heat balance equation 


can then be written 
Qep (T, T) t | 
U(F)(T 


, 


AH) kV AP 
T.) 


_oe 
dé 
(1.7) 


The controller varies the coolant flow rate in 
response to measured deviations of the tempera- 
ture from the desired value 7, and to the deriva- 
tive of the temperature variation. For reasons 
given earlier (Part II, section 9, [1]), 7’, is chosen 
to be a natural steady state of the system. Thus 
the controller equation can be written 


F = F,+¢,(T —T,) 4 cg (T - 7) 

It is this equation together with (1.1) and (1.7) 
which have been mtegrated to show the behaviour 
of the system, the interest again being in whether 
the naturally unstable steady state can be con- 
trolled, 
which have been chosen to simulate a realistic 


The values of the physical constants, 


situation, are given in Table 1. 


Physical constants 





Constant Value t nits 
1-000 Ib moles ft* 
650 RK 
os ft" se« 
60-0 Ib ft® 
62-5 Ib ft® 
10 B.t.u lb R 
0 B.t.u lb R 
100 re 


Constant Value Units 


500 ft? 
0 | B.tu ft? br oR 
100 | Btu ft? hr °R 
20,000 B.t.u. Ib mole 
27 » 10! ft? Ib mole sec 
520 *R 
4.700 B.t.u/lb mole 
11 B.t.u Ib mole “R 
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Table 2. 


Dimensionless parameters 








A steady state has concentration A = A, and 


temperature T’ — T, satisfying the equations 

Q(A, — A) — hk, VAZ = 0 

T,) AH) k, VAZ 
U,(T, 


(1.8) 


Qcp (T, 
T)=0 (1.9) 
where k, = ky exp(— E/RT,) and U, = U (F,). 
For the numerical values of Table 1 the unstable 
steady state we wish to control is at 

A, — 0-61912 Ib moles ft®, 7 
with «.. 


It is introduce 
quantities by the substitutions 


700-00 KR 


0-13810 ft® see. 


convenient to dimensionless 


tr=Q0V.a=(A-— A) A 


t=(T — T,)/T, f — F — F, 
With these substitutions equations (1,1) and (1.7) 
can be linearized in the usual way, and become 


da 


(l + aja xi yt (1.10) 


dr 


B—y)t— df (1.11) 
where the dimensionless constants are defined 
in Table 2 x, 8, y and 8 play the same role as in 
Part 1, section 3, though they differ in detail be- 
cause of the different definition of the dimension- 
less quantities. With this linearization we can 
discuss the stability of the steady states, 
2. Srapimiry Unper Conrro. 
The form of control we wish to investigate is 
that in which the departure of the coolant flow 


rate @. from its set point is proportional to a 
linear combination of the departure of the reactor 
temperature T from its set point 7, and the 
derivative of the temperature. In terms of the 
dimensionless co-ordinates introduced we put 


dt 


A 
dr 


Sf = M+, (2.1) 


where A and ,« are the two proportionality con- 
stants. Equations (1.10) and (1.11) then become 


da 
(l +a)a yt 


dr 
t (2.3) 


Proceeding exactly as in Part I, section 4 we have 
to examine the roots of the quadratic 


p+l1+z2 
A(p) = 


This may be written 
(1 + pe) p* +d,’ ' j 0 
where d,’ = (1+)(1 + ,) 


+ B—y+a (2.4) 


dy = (1+) (14 -A(1 +2) 


and the condition for stability is that (1 + ,), 
d, and d; should all be of the same sign. In the 
absence of control this means that 


+1+fB—y>0 and 
dy = (1 +2)(1 + B)—y>0 


the criteria obtained in Part I, section 3 equations 
(3.5) and (3.6) [1]. 
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We are interested however in controlling a 
naturally unstable state in which not both ¢, 


and d, are positive. The constants A and « there- 


fore have to make the three quantities 


d,+A+( a) pe, 


f | (1 a) A 


all of the same sign. It is intuitively evident that 
no good can come of large negative A and ,, for 
this would mean that the cooling rate would be 
decreased when the temperature increased, just 
the opposite of what is required. It would be a 
mistake however simply to require all three to be 
positive and dispose of the analytical possibility 
of all being negative with a flurry of intuitive 
arm-waving, for in the proper resolution of this 
difficulty we shall tind that not only ts y« 1 
(making all quantities positive) but also « 1, 
a restriction that could not otherwise have been 
anticipated. This will be the burden of the next 
paragraph ; we will conclude this one by seeing 
what features of the reactor’s behaviour can be 
predicted before the actual calculations are made. 

The values of the proportionality constants A 
and » may be taken as co-ordinates of a point in 
a plane. The region of stability defined by 
1+p>0,d/ > 0, d’ 
and to the right of the three lines AB. BC and 
CD in Fig. 1. 
are complex if the point A, «is within the parabola 
a," (1 


0 is thus the region above 
The roots of quadratic \(p) 0 


w)4,. This parabola touches AB 


(ita) Aeds+O 


SADDLE POINTS 


UNSTABLE NODES 


Fig. 1. Stability 


for control parameters 
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at B and CD at C and immediately shows the 
character of the steady state at any point, as 
shown in Fig. 1. Fig. 1 is drawn for the given 


system and the numerical results are 


A + 22300 » > 21695, A > 11855 


Since the steady states are detined by setting 
the derivative equal to zero it is evident that their 
position is independent of «. For a given value 
of A the position of all the critical points can be 
obtained by solving equations (1.8) and (1.9) for 
A, and T,. Moreover since A does not appear in 
equation (1.8) all critical points lie on the curve 
in the A, T 
The method of solution is exactly the same as in 
Part Il, section 9 [1] and need not be repeated 


plane defined by equation (1.8). 


here. The result of the calculation is shown in 
Fig. 2 where the numerical values are for the 
example given, As in the previous example it is 
found that the two other steady states coalesce 
and disappear at a value of A (1-46) less than that 
at which stability of the controlled steady state is 
assured. We may thus expect a region where there 
isa stable limit evele surrounding the still unstable 
control pont, 

\ more interesting situation however, which 
was uncovered in the previous work, may be 
expected here, namely the presence of an unstable 
limit eyele surrounding the control point even 
after this has been made stable. This means that 
the control is really only exercised over a small 


STABLE FOCI 


1 and “. 


[Numbered circles refer to figure number of phase portraits! 
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region of the phase plane and that any excitation 
that takes the system beyond this will lead to a 
steady oscillation round a stable limit cycle. 
This phenomenon is discussed in detail in Part IL, 
section 10 [1], but in the example there given 
was very transitory. In the present case we might 
expect a zone lying to the right and above the 
transition line BC such that values of A and «u 
lying there would give rise to this effect. Such 
is the case and the part of its boundary which has 
been determined is the line FGH in Fig. 1. 


CONCENTRATION 





TEMPERATURE 
oe 


600 700 





"s 


Fic, 2. Position of steady state for varying A. 


3. Srapsmrry wirt Detayep Conrrot 

If the controlling system has a small delay 
« in transmitting the signal the values of 7 and 
dT dé actuating the value for the coolant flow 
rate Q. at time @ will be T'(@ — «) and (dT d@),_,. 
Thus the transform of df will not now be 


if = M t(0)] but 


+ ph [ pi 
+ wp [pt — t(0)]} 


if exp ( (3.1) 


If this is inserted in the Laplace transformation 


of equations (1.10) and (1.11) and these are 
solved for / we have, 


ep) {At 


i= {«aa(0) +(p +1 +2) 
{1 + p exp (— ep))t(0)}/ A’ (p) 
where 


A’ (p) + exp (— ep) 


x) (mp + A) 


(? j d, p- dy) 


(p +1 (3.2) 


The stability will now depend on the sign of the 
real part of the zeros of A’ (p), the system being 
stable if these are all negative. 





Lip) - PLANE 
Ariv2 pre 
(STABLE) 


L(p) - PLANE 
Aviv2 pr? 
(UNSTABLE) 


L(—) - PLANE 
Ast peels? 
(UNSTABLE | 


Fic. 3. Nyquist diagrams. 


The existence of roots with positive real part 
may be shown by the Nyquist diagram. Consider 
the function 


p? +dp +d 


3.3 
(p +1 + «)(ep + A) (3.3) 


L(p) = exp(— ep) + 
This will map the p contour of Fig. 3 into a curve 
in the L plane such as is shown in the other 
diagrams of this figure. The imaginary axis 
p = tw maps on the curve AA’ B’ B which goes 
from the point A on the real axis of the L plane 
and goes rapidly into the circle L = exp (— tew) 
+ » ‘ for when w, and so |p|, is large the quotient 
of quadratics is close to »'. The point B on the 
imaginary axis of p will correspond to some point 
B on the circle of argument w,. Now as the 
imaginary part of p is held constant and Re (p) 
increases the quotient of quadratics will remain 
close to »' since |p| is large and arg exp (— «p) 
will remain constant (= w,) but the modulus of 
exp (— ep) will tend to zero, The image of BC in 
the p-plane is thus the radius BC of the circle in 
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L plane. As p traverses CDE the image point : ; 
delay « then the function whose zeros must be 
will describe an arbitrarily small curve about the ; 


+, ; examined ts 
point L = « *. The rest of the curve is a mirror 


mage of the earlier part in the real axis of the 1 pad p+ dy 


L plane. lp) < ap (— @) A (p +1 + «) G@(p) 


Now the value of ¢ affects the approach to the 

PHASE PORTRAIT 
Aeiass 

but however small « may be the circle will have wr 

/ @ wer cree 


circle and the number of times it is traversed, 


unit radius and centre « '. Thus whatever the 
behaviour of the rest of the curve the circle will 
certainly contain the origin if — 1 < » ' 1 ie. 
je 1. We have then a further condition for 
stability t< pl. The first mequality 
disposes of the paradoxical region ! 
A< d, (1 + 2), the second is a new restriction 
not anticipated by simple linearization. 

The third and fourth sections of Fig. 3 show 
two cases of instability. In (ii) it is due to the 
delay term, it bemg the circle that encloses the oO, 0, Phase portrait. 
origin; m (iv) the conditions are unstable even 


PHASE PORTRAIT 


without delay. The curves are drawn for the Tips 


constants of the given svstem and for « ol, sa pe ve 
. A STABLE LmwiT CYCLE 
but would searcely be different for any other fi rr @ warene Lest crcut 
suiliciently small value of « 0. 
it will be observed that a restriction of this 
kind will apply in the general case discussed in 


Part I, section 5 [1] where 


if A G(p)t 


K™ 


\ 


\ ‘ 
\ \ . \ 
\ 
\ \ 


and G(p) M(p) N(p) is a rational function 
of degree m in the numerator and n in the deno- ~ ~ 


minator, if m — nn + 1. Physically this transfer 
. » > » Ne @ 7 © > Phase . 
function defines an imperfect form of derivative Fig. G. A = 16.6 » = 0-25 Phase portrait. 
and proportional control. If in addition it has a 
PHASE PORTRAIT 
PHASE PORTRA! hens 
Aso we 
ONLY THE SEPARATRICES 
ARE SHOWN 


\ 


a a 2 a. 
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0,» ~ 0 Phase portrait. Fe 1-3, « = O35 Phase portrait. 
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Since m " 1 the second term is the quotient 


2 and sO 


_ pe 


of two polynomials each of degree 1 
exp ( ep) 
Again we must require |A) < 1, 


when |p| is very large L (p) 


t. Resuvrs 

As in the earlier parts of the paper the behaviour 
of the reactor has been studied by integrating the 
non-linear equations (1.1) and (1.7) numerically, 
The Univae Scientific Computer, Model 1103, 
was again used and the Gill modification of the 
Runge-Kutta numerical integration scheme. The 
general structure of the programme has been de- 
scribed in Part IIL.12 [1] and need not be repeated. 
The value of the derivative for control purposes 
was taken to be the value at the previous time 
increment. Thus the slight delay « was simulated 
and the 
|e > | 
effect on the calculation of curves in the phase 


regarding instability for 


That this delay had no 


conclusion 
was checked. 


plane was established by decreasing it to 1 32 of 


its value; with this decrease the curves were 
reproduced accurately. 

Of the phase portraits calculated Figs. 4-8 
will suflice to show the general behaviour of the 
system. Their the 
boundaries of stability is shown in Fig. 1. Fig. 4 
is the natural face of the reactor, showing the way 


This is entirely 


disposition relative to 


without control. 
similar to the system 
(Part II, Fig. 11). The naturally unstable steady 


state at T — 700 R, A 0-61912 moles cu. ft is 


it behaves 


previously discussed 


In Fig. 5 A 
control to 


the one we wish to control. 1-66, 
he 0 


removed the other states, but the control state 


and there is. suflicient have 
is still unstable and surrounded by a stable limit 
Fig. 6 (A 1-66, 0-25) the 


condition of partial control in which the stable 


evele, shows 


control point is surrounded by an unstable limit 
The 


unstable limit cycle was accurately located by 


evele and this again by a stable cycle. 


reversing the direction of time. This reverses 
the stabilities and by starting a trajectory very 
close to the control point, which ts now unstable, 
the solution spirals outwards to the now stable 
limit evele. By a careful search using this 
technique the boundary FGH (see Fig. 1) of the 
region of the A, « plane in which such unstable 
limit cycles exist was determined, The sharp 
vertical boundary FG seems to correspond to the 
value of A above which only one steady exists. 
This is confirmed by Fig. 7 with A = 13 y = 0-5; 
although this is very close to the threshold of 
stability there is no limit cycle but two unstable 
For clarity only the 
Finally 


0 shows the whole 


and one stable critical point. 
separatrices are shown in this figure. 
Fig. 8 with A = 3-33 and yu 

phase plane under full control. 


NoTATION 
concentration of reactant 
inlet concentration 
steady state concentration 
area of cooling coil 
activation energy 
2Q. Co ht Ay, 
steady state value of F 
M ip) N(p) 
control 


dimensionless flow rate 
transfer function for imperfect 


heat of reaction 

proportionality constant 

complex variable of Nyquist plots 
numerator of G (p) 

denominator of G (p) 

flow rate through reactor 

flow rate of coolant 

temperature 

inlet coolant temperature 

inlet temperature 

steady state temperature 

heat transfer coeflicient 

steady state heat transfer coeflicient 
volume of reactor 

(A Ay) A 


s 
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specific heat of reactant 

specific heat of coolant 
biadt4+B-y 

(1 x)(1 + B)- y 

d, (l + apsa 

dy a x)A 

F — F, 

heat transfer coeflicients for cooling coil 
reaction rate constant = k,exp(-— ER 


steady state value of k 


REFERENCE 


Laplace transform variable 
(T— T,) T, 

time delay in system 

time 

proportional control constant 
derivative control constant 
density of reactants 

density of coolant 

gov 


imaginary part of p 


{1} Anis R. and Amunpson N. R. Chem. Engng. Sci. 1958 7 121. 
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The hydrodynamics of flow between horizontal concentric cylinders 


Generalized criteria and correlation of heat transfer data 


A. H. and F. C. Haas 


Department of Chemical Engineering, Rensselaer Polytechnic Institute, Troy, New York 


NISSAN 


(Received 16 May 1959) 
Abstract The neutral stability criteria for the generalized case of flow due to circumferential 
pumping and rotation of the inner cylinder are reviewed. The neutral stability is defined by the 
parameter, P, as used by Brewster ef al. Heat transfer data taken from the literature have 
been recalculated. The parameter P is found to be useful in correlating the data. 


Résumé Les auteurs ont examiné les critéres de stabilité indifférente pour le cas généralisé 
dun écoulement da a un pompage circulaire et a une rotation d'un cylindre intérieur. La stabilité 
indifférente est détinie par le paramétre P, tel qu'il a été utilisé par Brewster ef al. Les données 
du transfert thermique prises dans la littérature ont été calculécs a nouveau. Hl a été constaté 
que le paramétre P a son utilité dans la corrélation des données, 


Zusammenfassung —-Verfasser teilen das neutrale Stabilitatskriterium mit fiir eine Strémung, 
die unter dem Eintluss einer Zwangskraft in Umfangsrichtung bei Retaticn des inneren Zylinders 


zustande komunt. 
detiniert. 


Wiarmeitibertragungswerte der Literatur werden von neuem berechnet. 


Die neutrale Stabilitét wird durch den Parameter P nach Brewster et al. 


Hierbei erweist 


sich der Parameter P als zweckmiissig, um die Daten zu korrelieren. 


INTRODUCTION 
Ix an earlier paper Brewster and Nissan [1] 
have discussed the fluid motion between concentric 


cylinders when the flow was due to rotation of 


the inner cylinder. The secondary flow patterns 
discussed included those that occurred when the 
primary flow was either fully or partially reversed. 
It is the purpose of this paper to review the case 
where the flow may be due to a circumferential 
pressure gradient and/or rotation of the inner 
cylinder and to show that a generalized parameter 
may be used to correlate heat transfer data. 
Review or GENERALIZED CRITERIA 

Dean [2] theoretically and 
Brewster et al. [3] have shown experimentally 
that when a fluid is pumped circumferentially 
around an annulus, a vortex instability occurs in 
the outer section of the annulus. The instability 
is similar to the vortex motion predicted and 
shown by TayLor [4] in his classical paper on 
hydrodynamics of fluids contained between con- 
centric rotating cylinders. 


has shown 


The problem of simultaneous circumferential 


pumping and rotation of the inner cylinder has 
been treated theoretically and experimentally by 
Nissan et al. [3] and theoretically by Di Prima 
[5]. It is observed that all possible combinations 


of secondary flow can be obtained by suitably 
adjusting the variables in this generalized case. 
To allow for the combined rotation and pumping, 

the parameter N is detined : 
na Mad (1) 

2q 
It may be noted from the definition of N that a 
value of zero corresponds to the problem Dran 
[2] considered and that a value of one corres- 
ponds to the problem considered by Tayor [4]. 
Nissan et al. [3] have considered the case of 
N | oo. Applying the Rayleigh instability 
criteria for inviscid fluids (shown by Synce [6] 
also to be applicable to viscous fluids) to this 
case, it is found that actually two regions of 
potential instability occur in the annulus. One 
region occurs adjacent to the inner cylinder wall, 
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to the 
evlinder wall. It is further shown that the region 
to the 
unstable first and this 


while the other occurs adjacent outer 
become 
the 
The 
criteria for neutral stability for the three cases 
N-O.1, 

If the 
streamline of zero velocity, the width is easily 


will 
to 


adjacent inner evlinder 


is considered he 
point of neutral stability for the system. 
and 4 © are listed in Table 1, 


unstable region, 5, is bounded by a 
caleulated by considering the radius of the zero 
velocity streamline as derived from the Couette 


flow equation ; 


fe \h a)*) 


92, (1 pe) 


1 (a b)? 


(2) 


a\* 22, (1 
("] 


-0, then 6 d may be detined as follows : 


(bay? 


If ad a 


Nissan ef al. (3) define a parameter, P?. which 
for all of N 


applied to the case of simultaneous circumferential 


deseribes the stability values as 


pumping and rotation of the inner cylinder. 


re) 


vy Ala 


(33) 


The values of P which detine the condition of 


neutral stability are given as: 


Oo< bd 23 12 


4 3 o d YO-7 (o dy> : 


10 
y [2(8 d) — 1] 

t 

y {[(2(6 d)— 1}? + O-01; (4) 


Table 1. Criteria 


Hear Transfer Data 


amount of data is 


for 


Although only a small 
in the literature 
through a fluid 


rotating cylinders, it is evident that the heat 


available heat transfer 


contained between concentric 
transfer rate is dependent upon the annulus size 
and the ratio of angular velocities of the cylinders, 
BaokkLUND Kays [7] the only 
extensive set of data at the present time. Their 
data indicate that the Nusselt number is constant 


and present 


and dependent only on geometry until the critical 
When the critical 
angular velocity is exceeded, the Nusselt number 


angular velocity is reached, 


steadily increases. The data indicate, as is to be 
expected, that each combination of da and 4 
presents an individual critical Taylor number. 
The authors use the following empirical equation 


to correlate their data : 
Nu 
Nuleond 


{Ta (Tae, 
| 41-4 + [Ta,, 


Tac, —o] [1 3.5 (d a)}\' . 


“ eta ol [3-5 (d a)| | 


(5) 


Gaz.iey [8] considers the problem of heat transfer 
in electric motors and presents a limited amount of 
heat transfer data for smooth concentric rotating 
evlinders. 


CORRELATION OF HEAT TRANSFER DATA 
The data of Baorkiuunn and Kays [7] and 
Becker [9] have been recalculated using the 


parameter P. Since there was no circumferential 


for neutral stability 





Defining equation 


Width of unstable region | Author 


| Dran [2] 


| Tayior [4] 


Nissan ef al. (3) 
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pumping involved, the parameter may be detined 
as: 
2,a5 


P (6) 


Vs a) 
or more conveniently for this case : 


P 


1/2(5 d)'/* Ta, 


where Ta, is the Taylor number which has been 
corrected for the finite gap width by the use of a 
geometric factor, F,, as suggested by Kaye and 
EiGar [10]. The geometric factor modifies the 
Taylor number when the gap width is large, ie., 
when d a does not approach zero. 


Ta = Ta, F 


z£ 


where F, is detined by 


ad Tm 


br.) 


i?» 
b -! m 


0-0056 | l 0-652 (, (b lm) 


)| 


Critical values of P as 


Oy, 
b -! m 


Table 


For counter-rotation of the two cylinders, 6 
calculated radii of 
streamline velocity through use of the Couette 


was by considering zero 


flow equation, 
b?) 
-p(b a)*] 


(ju 1) (a? 
(a by?) [1 


For co-rotation of the two cylinders the unstable 
region extends over the entire gap width. Follow- 
ing the suggestion of DiPrima [11] the angular 
velocity was approximated by a constant, the 
arithmetic mean velocity, and the parameter was 


detined as: om 
ra (1 


' 


P #) 

The points of neutral stability defined by P. 
as determined from the data are listed in Table 2. 
23 


1 are quite constant and agree well 


It is noteworthy that the values for 6d .- 
and 6 d 
with the theoretical values previously cited. The 
data have been plotted using the parameter 
P P. to represent the hydrodynamic situation 
(Fig. 1). The following equations may be said to 
represent the data. 


determined from heat transfer data 





da 


0-054 
O-OS4 
O-Os4 
-Os4 
0-128 
(128 
128 
O24 
Ol 
2 


ooo 
O56 
ooo 
16 
O55 
O06 
1-06 
O55 
O06 


105 


O54 
O-OS84 
0-084 
128 
0-128 
O24 
O24 


Reference 


oh ef ef oS eh aS aS oS of aS 





— 
~ 
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Fie. 1 


d Ry 


: 


30 40 6©HO 


Correlation of heat transfer data using the parameter PP, 


Author 


Brorktunp and Kays 
Boorktunp and Kays 
Beorkiunp and Kays 
BsoRKLUND and Kays 


246 
0-128 
O-O84 
moss 


Lae 3) 


Nu P 
N Uoond | P. 


os 


0.5 


: Nu P 
3 « =< 40 ; (p 


; 
Nukond. 


All experimental points lie within 15 per cent 
of the lines detined by the above equations. 
Within the limits of present knowledge (all data 
available is for one tluid-air) the above correlation 
allows one to predict the heat transfer charae- 
teristics of a concentric rotating 
cylinders if the geometry and cylinder rotational 


The value of P. is given 


system of 


speeds are specified. 
directly by a knowledge of 5 which may be 
evaluated from the geometry and «. P is evaluated 
from knowledge of the cylinder rotational speed, 
and with these the heat 
transfer coefficient may be evaluated. 

It would be desirable to know the effects of 


values, convection 


Becker 


Prandt! number and circumferential pumping on 
the are not 
available at present. 


above correlation, but such data 


CONCLUSION 


It has been established that the parameter P 
is useful in correlating heat transfer data for a 
fluid contained rotating 
evlinders. A study of data taken from the 
literature has indicated that the convection heat 
transfer the 
geometry and cylinder angular velocities are 


between concentric 


coctlicient may be predicted if 


established. 
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NOTATION 


radius of inner cylinder radius of zero velocity streamlines 
radius of outer cylinder axial length of annular space 
width of annulus "= average velocity in unstable region 
geometric factor detined by equation (7) 6 = width of unstable flow region 
convection heat transfer coeflicient ratio of angular velocity of outer and inner 
thermal conductivity cylinders 

= dimensionless number defined by equation (1) » = kinematic viscosity 

= 2hd k = Nusselt Number angular velocity 

2(d a) 

In[(d/a) + 1] 
duction 
(2, ad)/v y/(d/a) = Taylor number 
dimensionless number defined by equation (3) 
volumetric flow rate through pump arithmetic mean 
radius denotes cylinder speed ratio 


~ 2 Nusselt number for con- Subscripts 
denotes d/a +0 
denotes inner cylinder 


denotes condition of neutral stability 
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Heat transfer in two-phase flow 
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Abstract Heat transfer measurements for the twoephose flow of water air and gas-oil air 


muatures in a vertical tube are deseribed 


Much attention has been given to the contribution of the heat-transfer section of the tube 
and the thermocouples for measurements of tts wall temperatures so as to obtain accurate results, 
re compared with those for pure liquid How. Dt is shown 


The results of the measurements a 


that especially the first amount of gas added to the liquid raises the heat transfer coeflicient 


considerably, The results for the two-phase systems can be represented by simple correlations 


between the Nusselt and Revoolds numbers of the svstens 


\t higher gas liquid rotios a maximum im the heat transfer coeflictent is observed. An 


attempt is made to characterize its position by a dimensionless quantity 


Résumeée On a frit des mesures du transfert de la chlewr dans le cas d'ecoulement de svstemes 


a deux phases, a savoir eau air et gas-oil air, dans un tube vertical 


Pour obtenir des résultats exacts. une grande attention a ete devouee a la construction 


dle la section du tube ou le transfert de chaleur se deroule, et aux thermocouples servant a mesurer 


les temperatures a la parot 


Les resultats obtenus ont ete compares avee ceux pour Tecoulement dun liquide normal. 
Ila ete demontre que c'est surtout la premiere quantite de gaz ajoutee au liquide qui donne lieu 


a une forte augmentation du coeflicrent de transfert de la chaleur, Les resultats pour les systemes 


a deux phases peuvent étre representes par de simples correlations entre les nombres de Nusselt 


et de Revonolds des svstéemes 


\ des taux gaz liquide plus cleves on observe un maximum dans le coeflicient de transfert 


de chaleur, On a fait un effort de carectertser la position de ce maximum a Taide dune quantite 


seers CLINE TSE oth 


Zusammentassung) Dice Verfosscr beschreiben Wiarmetibergangsmessungen fur die Zwei- 


phasenstromung von Wosser Luft and Gas-OL Luft in einem senkreehten Rober 


Auf die Konstruktion der Messstrecke des Rolres und der Thermoelemente zur Messung 


der Wandtemperatur wurde besondere Sorgfalt’ verwendet 


Die Messerygebnisse werden mut denen fair reine Flissigkeitsstromung verglichen. Es zeigt 
Gasmenge, die der Flissigkeit zugefiigt wird, den Warmeitiber 


sich, dass besonders die erste 
Die Ergebnisse fir das Zweiphosensystem konnen durch 


gangskoeflizienten betrachtlich erhoht 
eine cinfache Bezichung eawischen den Nusselt- und den Reyvnolds-Zahlen des Systems dargestellt 


werden 


Bei hoheren Verhiltnissen Gas Plissigkeit wird cin Maximum des Wiarmetibergangs- 


koeflizienten beobachtet Seine Lawe wird versuchsweise durch cine dimensionslose Grosse 


charakterisiert 





Heat transfer in two-phase flow 


1. INrrRopUcTION 


Ix 1951, Sremerpinc and Vexrscnoor 
published results of measurements on 

transfer to water air mixtures 
currently through a vertical tube. Their imme- 


diate purpose was to obtain data for the simul- 


(1) 
heat 


flowing con- 


taneous heating of a liquid and a permanent gas, 
to be used for the design of a heat exchanger 
for heating up a liquid gas mixture to reaction 
temperature. They the liquid 
velocity and the air-to-liquid ratio and found 
that over a wide range the coeflicient of heat 


varied mass- 


transfer increases with increasing values of these 
quantities, 

After the publication of their results on the 
system water air, the work in this laboratory 
was continued with measurements on gas-oil ‘air 
mixtures. It then appeared that in some cases 
the measurements of wall temperatures were not 
sufliciently accurate. An analysis of this difficulty 
resulted in the construction of a new measuring 
tube of considerably increased wall thickness 
and with improved mounting of the thermocouples. 

In this article the considerations which led to 
the improved construction are presented as well 
as the results obtained, both with water air and 
gas-oil air mixtures. In a general way these 
confirm the trends observed by the previous 
authors. 

2. DESIGN OF APPARATUS 
The heat transfer coeflicient from a_ heated 
tube wall to a liquid flowing inside is given by the 
well-known equation : 


te = % (ty — ty) (1) 


¢, = heat transferred per unit time and unit 
area, 
x = heat transfer coetlicient 
t,, = wall temperature 
t, = bulk temperature of liquid. 


Of the variables in this equation only ¢,, the 
temperature of the wall where it is in contact 
with the fluid to be heated, is difficult to measure. 

For a correct determination of the average wall 
temperature two conditions must be fulfilled. 
In the first place the insertion of the sensing 


element, in this case a thermocouple, in the wall 
must be such that the heat flow through the 
wall is not perceptibly disturbed. Otherwise the 
temperature at the measuring point can adjust 
itself to a value differing from that in undisturbed 
surroundings. In a thin-walled tube this condition 
fulfilled. Secondly, the 
temperature distribution along the wall should be 
as uniform as possible. However, when saturated 


cannot practically be 


steam is used as a heating medium, considerable 
local temperature differences occur on the heated 
side. For, owing to differences in wettability 
of the wall the condensate runs down in bands 
and the rest of the wall remains dry. Thus the 
dry and the wetted spots will develop temperature 
differences which are of the same order of magni- 
tude as the temperature drop across the condensate 
film. Under certain conditions the latter was 
sometimes 20 C in our experiments. In a thin- 
walled tube these differences are only partly 
levelled out by conduction and the result is an 
irregular temperature distribution on the side 
where the coetlicient of heat transfer is measured. 
Consequently the average temperature on that 
side cannot be determined accurately. 


“ap 
Heat remow ai 








VAAN 


a 
Y 


> 


x 


Fic. 1. Heat transport through a thick wall. 


In order to obtain an impression of the wall- 
thickness required to smooth out these tempera- 
ture differences sufficiently, a calculation was 
made for the following simplified model: (see 
Fig. 1). 

An infinitely long wall of thickness d is heated 
on the side x = 0. On that side the temperature 
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shows a sinusoidal disturbance in the y-direction 
(parallel to the wall), so: 
Y 
(0, y) = t, + Atcos 27 7 (2 

At is the amplitude of the disturbance, / the 
distance between two temperature maxima. 

Ata d the heat is transferred in accordance 
with equation (1), so that: 


(3) 


y id) i 
hdd x [t(d, y) — ty | 


A stands for the coetlicient of heat conductivity 

of the wall material. For the stationary case we 
have: 

ot . Mt 

’ (4) 


” ” 
d rr y- 


Substituting a d in the solution of this set of 
equations the following expression is found for the 
temperature distribution at the side of the wall 
d): 


where heat is removed (a 


(tf t,)axd A 
2 (d, y) : 
\ ’ iL (add \) | 


Qnrdsl = 
(2rd l)sinh2ad } 


tied \) 
Atcos 27 yl 


t N 


° R At cos 27 y / 

In this equation N represents the temperature 
drop across the wall in the absence of a distur- 
bance and R denotes the fraction of the distur- 
bance amplitude Af remaining on the opposite 
In Fig. 2 R has been plotted 


two 


side of the wall. 
ld 


to wall 


between 
A ad 


resistance in liquid film to resistance in tube wall) 


(ratio of distance 


thickness) with 


against 
maxima (ratio of 
as parameter, 

It is seen from this Fig. that in order to reduce 
a disturbance of 20 C at the heated side to a 
value of 2. C at the liquid side of the tube, a 
value of 2 for ld is the allowed. 
The distance / is diflicult to estimate, but for a 


tube of 14mm internal diameter, 10 mm seemed 


maximum 


reasonable in view of visual observations in a 


similar, transparent tube and this leads to a wall 


thickness of 5mm. The material of the wall is 


HENDAL 


not important in this respect, since A/ad has 
little influence on R. The choice fell on copperin 
order to keep down the temperature gradient in 
the radial direction. In with the 
location of thermocouples the wall was made 
slightly thicker, namely 6 mm. 


connexion 


| 


| 


A - Resistance in liquid 


of termpercture disturbance 


oo Resstonce in tube wo! 





Remaining fraction 


moe between 


Wall theckr 


temperature mon 


Fic. 2. Smoothing-out of temperature disturbances 


in wall of heat transfer tube. 


A correct insertion of the thermocouples is 
facilitated by this wall 
Two constructions were used, which 


greatly considerable 
thickness. 


are shown in Figs. 3(a) and 3(b). 


In case (a) a hole of 2mm diameter was drilled in 


the 
\ 0235 mm constantan wire was attached to a thin 


wall, reaching to 1-3mm from the inner wall. 


conical copper disk 2mm in diameter, This was 
placed in the hole and subsequently pressed flat, thus 
ensuring a good contact with the bottom. The rest 
of the hole was then filled with close-fitting copper 
rings with a bore of 1mm, leaving space to insulate 
the wire with a porcelain tube. In (b) a 1mm hole 
was drilled in the wall and the porcelain-insulated 
constantan wire was introduced into it. Perpendicular 
to this hole a second one was drilled, into which a 
close-fitting copper evlinder was driven to fix the wire. 


214 





Heat transfer in two-phase flow 


Fic. 3. Methods of attaching thermocouples in tube 


wall. 


In both the wire forms, 


together with the tube material, a copper con- 


Cases constantan 
stantan thermocouple with a well-defined junction, 
which does not perceptibly interfere with the heat 
low through the wall. Four couples of each type 
were installed along the tube wall, two of different 
The 


couples were calibrated at two fixed temperatures. 


types at the same height in each case. 


During the experiments the readings of the 
couples at the same height differed by not more 
than 0.5 °C, while the temperature profile in the 


axial direction was very regular. It may, there- 


fore, be assumed that with this combination of 


the 
average wall temperature is measured with an 


wall thickness and thermocouple design 


accuracy better than 1 °C, 


The measuring tube itself had a length of 


20 cm, internal diameter 1-4em and was placed 
in a double-walled steam jacket, as shown in 
Fig. 4. Steam was introduced at three places 
(A) in the outer jacket and reached the measuring 
section through an annular opening at the bottom 
of the inner jacket. The condensate formed on 
the measuring tube was removed through B and 
its quantity and temperature were measured, 
Condensate on the outer jacket was removed 
through C, whilst the amount of steam supplied 
was so adjusted that a small excess escaped at D. 


Excess steam 
4 
| 

8) 


Thermal insuia tion 
(Epikote’) 
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Copper measuring tube 
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Condensate from 
outer jocket 


Condensate from 


easurng tube Non-neoted section 


ength :i20 cm 








Directon of flow 


Apparatus for measuring coeflicients of heat 
transfer. 


the axial 


would be 


In order to avoid heat losses in 
the 


siderable at the wall-thickness applied, Epikote 


direction of tube, which con- 
synthetic resin plates of 22 mm thickness provided 
with a hole exactly 14mm in diameter were 
placed between the measuring tube and the inlet 
and outlet tubes. The non-heated section of the 
tube was 120 cm long. 

The following quantities were measured ; tem- 
perature of liquid entering and leaving the tube 
(after separation of gas), the temperature and 
volume of the condensate formed during a certain 
period and the wall temperature at 8 points. 
The difference in wall temperature between top 
and bottom of the measuring tube never exceeded 
8 °C; as this difference is small with respect to 
(t,, — 4), the arithmetical mean of these 8 tem- 
peratures could be used in the calculations. By 
parallel connexion of the thermocouples, all of 


which had the same internal resistance, this mean 
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was read from a Brown self-balancing potentio- 
meter. It was then corrected for the temperature 
gradient in the radial direction in the tube wall, 
calculated from the total heat flow, the surface 
area and the thermal conductivity of the material. 
This correction never exceeded 1.5 “C, The values 
of « were calculated with formula (1), in which 4, 
could be taken as the arithmetical mean of liquid 
inlet and outlet temperature. 

ron Warer 


Air 


The apparatus was tested by carrying out 


Resvutrs or MEASUREMENTS 


WD Gas-O1 Wreriuout 


measurements on water and gas-oil and comparing 
the results with relationships known from litera- 
ture. In order to prevent the formation of air 
the the 


deaerated, and recycled in a practically closed 


bubbles on tube wall liquids were 
system. This included a cooler and was connected 
with the atmosphere through a narrow, liquid- 
filled tube. The deposition of scale on the heated 
section was prevented by the use of distilled 
water, while 0-2°, wt of potassium bichromate 
was added to prevent corrosion in the iron part 
of the test installation. 
The results of the measurements are given in 
Fig. 5, where 
; Nu 
"6 Pr'* \% 
eee eee 


"ee y" i 


‘ 
+? 
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Fic. 5. Correlation of heat transfer data for water 
and gas-oil in single-phase flow. 


Henpat 


Here » is the 
viscosity of the bulk of the liquid, », the viscosity 
For Re > 5000, the 
experimental points for water lic on a straight 


has been plotted against log Re. 
at the wall temperature. 


line, representing the relation : 
Nu = 0.030 Re? ™ Pr'/® (yn, /y,)° 


4000) lie 
on a straight line having the same slope, but the 


For gas-oil the points (provided Re 


coetlicient in this case is 0-028. This difference has 
little significance, because there is a possible error 
of 10°, in the values for the thermal conductivity 
and the specific heat of the gas oil, It is seen from 
Fig. 5 that the experimental values of the Nu 
numbers for water are about 20°, higher than 
predicted by the equation of Simprn and Tare [2]; 


Nu 


0-027 Re®*® Pr' * (m, »,,)""* 


This can be partly explained by the relatively 
small length-to-diameter ratio of the test-section, 
Theoretical considerations by Detssuer [3] which 
were experimentally contirmed by Harrnerr [4] 
have shown that over approximately the first 
five diameters’ length of a heat-transfer tube the 
thermal boundary layer is not yet fully developed 
and this results in high values for the heat transfer 
coetlicients, For our experimental conditions an 
increase of 10°, above the asymptotic value of 
the Nu numbers obtained in a long tube is in 
agreement with these considerations. 

The repeatability of the measurements 
amounted = to 3°... The of the 
experimental data is smaller than that of most 
literature data. From this fact and from the 
correlation observed it was concluded that the 


about scatter 


results obtained with this apparatus are reliable. 


Ain, Warer 
MIXTURES 


MEASUREMENTS ON 
Air, Gas-O1L 


(a) Air /water 


AND 


Air and water were supplied separately through 
needle at the of the 
unheated section of the tube at angles of 45° to 
To eliminate the effect of pressure 
pulsations in the tube on pump capacity and 
flow measuring instruments, a gear pump was 
used for the water supply and a large pressure 


valves mounted bottom 


its axis. 
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drop was allowed across the valves. With these 
provisions the flows became very regular. 

The mass velocity M of the water was varied 
from about 20 to about the 
volumetric air, water-ratio (V, Vj) from one to 


80 g cm? sec™!, 


beyond the value where the coetlicient of heat 


transfer reaches a maximum. Because it was 
found that the wall temperature influenced the 
results 
the 
coetlicients of heat transfer have been recalculated 
for a wall temperature of 60°C, The pressure 


drop over the water inlet valve was not kept 


5°, for a temperature rise of 10°C all 


constant in these experiments, but no differences 
were found between duplicate measurements. 
The repeatability was better than 10°,. The 
results of the experiments have been plotted in 
Fig. 6. 

It appears from Fig. 6 that for each value of M 
there is an unmistakable maximum in the value 
of x. As M increases, the value of V, J) at 
which the maximum occurs decreases. 




















6 8100 


Fic. 6. Heat transfer measurements on water air 


mixtures. 


(b) Air gas-oil 

Fig. 7 gives the results of the measurements on 
this system for mass velocities of 22, 31 and 
42 g em~? sec”! (properties of the gas-oil : 


0-83 g em~4 

23 dyn em! at 97 C 
2-40 ¢P at 15 C 
0-84 cP at 97 °C). 


The shape of the « versus V/V, curves is almost 
equal to that for water air, Again for each value 
of M a maximum value of « is found, though at 


the value of « rising by an average of 


considerable lower air liquid ratios than with the 
system water/air. 

















hic. 7. Heat transfer measurements on gas-oil ‘air 


mixtures, 


In this case it appeared that the value of « 
is influenced by the pressure drop over the gas- 
oil inlet valve, especially at low gas ‘liquid ratios. 
Raising the initial pressure from 1 to 2 atg. caused 
an increase in « values of 24°, for the experiments 
with VV, = 1, of 10°, when V, 
the effect was negligible for a gas /liquid ratio 
of 20. 
due to the better and finer mixing of gas and 
liquid at the higher pressure. A V,/V, = 20 a 
very fine dispersion is already obtained, caused 


V; = 5, while 


Photographs showed that the effect is 


by the mixing action of the air and the effect is 
consequently small at this ratio. In view of this 
effect the initial pressure was kept constant at 
1-6 atg in all gas-oil air experiments. The wall 
temperature also was constant at ~ 95 °C, making 
temperature corrections unnecessary. 

RESULTS 


5. CORRELATION OF EXPERIMENTAL 


To obtain a better impression of the influence 
of air addition on the coeflicient of heat transfer, 
the quotient «,/«, was plotted against the air-to- 
liquid ratio V,/V, for the conditions below the 
maximum in the heat transfer coefficient (Fig. 8). 
Here x, is the coeflicient of heat transfer for the 
two-phase systems, «, that for the liquids in 
single flow at the same values of the mass velocity 
M. 

For x, the experimental values have been used. 
Only for the gas-oil measurements at Re > 3000, «, 
was calculated by means of the relation deduced 
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coeflicient of heat transfer in 2-phase system 


os 
Ratio — = 
2, coefficient of heat transfer in single liquid flow 

(at the same mass velocity) as a function of air liquid 


ratho 


by Sreper and Tarr [2] for viscous substances in 
the range of laminar flow : 
(6) 


Nu 1-86 


d 
Re P 
¢ aly 


where d — diameter and /] = length of the tube. 

It is evident from Fig. 8 that the first amount of 
air causes a rapid increase in the coetlicient of 
heat transfer, but that for a further rise a very 
large quantity is required, Further it 


that the influence of air is most pronounced at the 


is clear 


lowest Re numbers. As the resistance agaist 
heat transfer is mainly due to the presence of a 
viscous sub-layer at the wall, this can be readily 
explained by the fact that the eddies, present 
in the wake of the rising air bubbles, penetrate 
into this viscous sub-layer and reduce its effective 
thickness. The thicker the original laminar layer, 
i.e. the lower the Re number, the more pronounced 
will be the effect which can be expected from this 
extra source of turbulence. 

For each of the systems used, a satisfactory 
correlation could be obtained by inter-relating 
the Nusselt, Prandtl and Reynolds numbers for 
the two-phase systems. In this correlation the 


Nu and Pr numbers were based on the physical 


HENDAL 


properties of the liquid, because the heat transfer 
will be determined by the properties of the liquid, 
For the 
Re number a volumetric mean Re, was used for 


at least as long as it covers the wall. 


the two-phase systems, defined by : 


y Pe D (1 — hy Pe v, D 


™" "e 


Re, (7) 
where / is the fractional hold-up of liquid, vy, and 
vr, the actual mean velocities of liquid and gas. 
As the actual velocity is equal to the superficial 
velocity (subscript: s) divided by the fractional 
hold-up, equation (7) can be written as: 


Pi Is D 


UT] 


Re, Re, + Re,,. (8) 
In consequence Re, is obtained by adding the 
liquid and gas Reynolds numbers, both based on 
superticial This definition of Re, 
fulfills the condition that for the extreme values 
of VL Vv, (0) and 
determined by respectively liquid and gas alone. 


velocities. 
©) the tlow properties are 
In Fig. 9 the data for the two-phase systems 


(for V, V, 


way as indicated in Fig. 5 for single-phase flow. 


1) have been plotted in the same 














Correlation of heat transfer data for the two- 
phase systems water air and gas oil air. 


Fig. 9. 
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Up to just below the maximum value of Nu all 
measurements for water, air satisfy to within 15°, 
the correlation : 


Nu = 0.029 Re," Pr'/* (m, ‘»,,)°"* (9) 


and for gas-oil air : 


Nu = 2-6 Re,°™ Pr'? (m,/4,.)""* (10) 


For all values of the liquid mass velocity 
Re,, covered a range from about zero to a 
multiple of Re,,. 

In the gas-oil /air measurements the Re number 
for the liquid alone (Re,,) varied from 1400 to 
3500, for water/air it was invariably above 5000, 
It is remarkable that for the two-phase measure- 
ments with gas-oil air the exponent of Re, is 
close to the value of 0-33 for laminar single-phase 
flow (equation 6), whilst for water air it is of 
the same order of magnitude as with turbulent 
flow of a liquid alone, This agreement of the 
exponent with that of single-phase flow, which is 
illustrated in Fig. 10, might indicate that the 
gas-oil in this case is still in a “ laminar” 
condition, 





©) 


SL. at 
aCe a ew 
~_+> > anand 


+ + + bee 





+ + + e+ 


10000 
Re, .+Re,, 














+ e 


Fig. transfer data 


for single- and two-phase flow. 


10. Comparison between heat 


THE COEFFICIENT Heat 


TRANSFER 


6. MAXIMA IN OF 


Description of the maximum in the coeflicient 
of heat transfer is difficult without exact informa- 


two-phase flow 


tion about the flow pattern, which could not 
in the An 
obvious explanation would be to connect this 
maximum with the transition of froth or slug 


be observed present apparatus. 


flow into mist-annular flow. Now it was observed 
that beyond the above maximum the outlet 
temperature of the gas was higher than that of 
the liquid. This experimental fact suggests that 
the maximum in the 
might coincide with the flow condition at which, 


heat transfer coetlicient 
owing to the highly turbulent motion of the gas- 
liquid mixture with increasing amount of gas, 
randomly disturbed “dry ~ spots appear at the 
walls. The heat transfer beyond the maximum 
will then no longer be controlled by the heat 
transfer from wall to turbulently flowing liquid 
alone, but also by the transfer from wall to gas 
at the dry spots. 

On the basis of this hypothesis a criterion was 
derived for the at the 
liquid film on the wall is partly and temporarily 


flow conditions which 
disintegrated by the shearing stresses exerted by 
the flowt. This was assumed to happen if the 
Weber number exceeds a certain value ; therefore : 


pg t,? 8 
Y 


= constant at the maxima 


5 represents the thickness of the viscous film and 
follows from : 


u® 6 


Pm 
" 


constant u* — the shearing velocity 
5 and u* can be eliminated from these equations 


by means of the relations : 
shearing stress at the wall 
Pm 
Strat SF 


friction factor 


C (Pm Um D/ny'* (Blasius equation) 
Index m indicates pro- 
perties of the gas liquid 
mixture 


+The authors are indebted to Professor J. O. Hinze for 
ths helpful suggestions. 
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If our assumption is correct, the maxima can 


the dimensionless 


thus be characterized by 


criterion : 
constant, 


M 


of gas and liquid 


Ps v,* " | M, a | 5 
M, Y " 


, ~ total mass velocity 


The numerical values of this dimensionless 
quantity at the maxima are listed in the Table 
below. 


The average wall temperature at the maxima 


HenpatL 


was 60 °C for the water air experiments, 97 C 
for the gas-oil air system, 

The result is fairly satisfactory but it should be 
kept in mind that it is obtained from observations 
on two liquids only and that a number of pro- 
perties, e.g. the density of the gas, were not 
varied, Moreover the influence of the heat flix 
on the stability of the wall-film was neglected. 
This result, therefore, should be considered with 
some reserve as long as no further experimental 
evidence is available to contirm the hypothesis 


underlying the criterion. 
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NOTATION 


Nusselt number 


(*) 


t 
Prandtl number | , 


\ 
pa °) 


Reynolds number | 


’ 


Reynolds numbers defined by equations (7) 
amd (8) 

specific heat 
wall thickness om 


tube diameter om 
1 


eal g ' deg! 


” 
mass Velocity g om * see 
total mass velocity of gas and liquid 


gem? sec! 


temperature 

linear velocity eM se4 

superticial velocity em sec 

volumetric gas liquid ratio 

coefficient of heat transfer 

cal em-? deg”! sec 
dynem 

1 


surface tension 


dynamic viscosity gem! see 
coefficient of heat conductivity 

cal em! deg! sec 
density gen 


Subscripts : 
g and / are used for gas and liquid respectively 
w and + are used for wall and bulk (temperature and 


Viscosity) respectively. 
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1958 Heat Transfer and Fluid Mechanics Institute. 
Preprints of Papers held at University of California, 
Berkeley, California, June 19, 20, 21, 1958. Stanford 
University Press, Standford and Oxford University Press, 


London. viii . 264 pp. £3 $8.50, 


Tue Heat Transfer and Fluid Mechanics Institute is “a 
nonprofit-making organization whose primary purposes 
are to promote the preparation of, and provide facilities 
for, the presentation of scientific and technical work of 
a fundamental nature in thermodynamics, fluid mechanics 
The Institute 


was formed in 1948 and has since met once yearly for a 


and related fields of scientific endeavour.” 


three-day Symposium at one of its sponsoring Californian 
Universities for the presentation of selected papers in the 
above fields. These papers are made available as bound 
collections of “ Preprints,” and judging by the steady 
increase in size of these volumes, the first decade of the 
Institute's existence has been eminently successful. 

The volume under review is the eleventh in the series 
and, printed on a somewhat reduced format, has been 
offered at a slightly lower price than in the previous year. 
Twenty-one papers are included; they cover a wide 
range of topics and are broadly classified under the headings 
of Fluid Dynamics, Dynamics of Reactive Fluids, Heat 
Transfer, Magneto-hydrodynamics and Hypersonics, plus 
two survey articles on Theoretical Acrothermodynamics 


and Magnetohydrodynamics, It is not possible to review 


Books 


Annual Reports of the Progress of Chemistry Vol. 55 1958, 


T. A. 


GrissMAN: Principles of Organic Chemistry. 


them in detail here but to make some general remarks 
instead. 

Attention is directed to the statements in the preface 
that the preprints should not be considered a formal 
publication and that “the preprints are made available 
only as a basis for discussion and not as contributions in 
final form to the permanent literature.” The Institute 
has, therefore, issued its warning, although it is difficult 
not to consider a book costing sixty shillings, and generally 
available, a formal publication. For those attending the 
Symposium the preprints obviously serve their purpose. 
To most others the volume will appear sometime later as a 
library copy and will be regarded as a work of reference 
containing papers dealing with the early development of a 
wide assortment of lines of research. 

If it were possible to do so, the publication would 
benefit from the inclusion of a couple of introductory 
review articles giving a general picture of recent develop- 
ments in heat transfer and fluid mechanics, in preference 
to the existing two surveys which are essentially statements 
of general principles in their respective fields. Preprints 
or otherwise, it is reasonable to expect all mathematical 
equations either to be typewritten or to be handwritten 
more neatly than appear for example on pp. 2 and 109, 
\ glance at pp. 46 and 32 in the 1956 and 1957 issues 
indicate an acceptable standard of 


respectively will 


handwriting. 


J. O. Witkes 


RECEIVED 


The Chemical Society, London 1959. 


Freeman, London and San Francisco 1959. 
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